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AUTOMATED  TRANSFER  FUNCTION  MEASUREMENTS 

1.  INTRODUCTION 

Frequency  response  or  transfer  function  measurements  are 
important  to  component  and  circuit  performance.  The 
frequency  response  of  a  linear  system  is  a  frequency 
dependent  relation,  in  both  gain  and  phase  difference, 
between  the  steady  state  sinusoidal  input  and  output  of  a 
circuit   or   device. 

Oscilloscopes  can  be  used  to  make  gain  and  phase-shift 
measurements  because  of  their  voltage  and  timing 
properties.  With  the  oscilloscope,  the  user  can  observe  the 
waveforms  under  test.  The  gain  of  the  circuit  is  determined 
by  visually  measuring  the  peak  amplitudes  of  the  input  and 
output  waveforms  and   calculated   using   the   relation 

Gain(dB)  =  20  log10 (Vout/Vin) 
The  phase-shift  between  the  two  equivalent  frequencies  is 
based  on  the  difference  between  the  zero  crossings  in  a 
particular  direction  of  the  two  signals  and  can  be  measured 
by  the  dual  trace  method  or  Lissajous  pattern  method. 
Measurements  made  on  a  conventional  oscilloscope  require 
repeated  setup,  user  interpretation  of  the  CRT  screen,  and 
calculations.  Such  measurements  consume  time  and  suffer 
from   low   repeatibility. 

It  is  therefore  necessary  to  automate  these  measurements 
for    more    consistency,     accuracy    and    speed     than    manual 
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methods.  Several  automated  techniques  are  discussed  in  the 
literature.  This  thesis  discusses  some  techniques  for 
automating  frequency  response  measurements  using  computer 
controlled  lock-in  systems.  Also,  a  software  simulation  of 
the  measurement  procedures  is  described  and  compared  with 
the  results  obtained  earlier. 

A  brief  introduction  to  lock-in  systems  is  presented  in 
section  2.  Section  3  deals  with  the  measurement  procedures, 
practical  implementation  of  the  procedures  on  two  test 
circuits  and  probable  error  sources  in  the  measurements. 
Section  4  deals  with  the  software  techniques  for  frequency 
response  measurements  based  on  the  input  and  output  samples 
of  a  test  circuit. 


2.  LOCK-IN  SYSTEMS 

Transfer  function  measurements  may  involve  the 
measurement  of  signals  obscured  by  high  levels  of  noise 
and  interference.  In  many  experiments,  the  noise  level  due 
to  thermal  noise  alone  may  be  several  orders  greater  than 
the  signal  of  interest.  The  amplitude  and  phase  variations 
introduced  by  the  circuit  or  device  have  to  be  estimated 
under   adverse   conditions   of   signal-to-noise   ratio. 

The  problem  is  essentially  one  of  'signal  recovery1. 
This  problem  can  be  better  understood  by  an  examination 
of  the  signal  and  noise  voltages  that  appear  at  the  output 
of  a  typical  experimental  system.  A  frequency  domain 
transformation  of  the  output  helps  in  identifying  the 
signal  of  interest,  noise  and  interference  which  are 
obscured  in  the  time  domain  as  shown  in  fig  1.  In 
evaluating  the  effect  of  noise  on  a  signal,  the 
distribution  of  noise  components  with  frequency  is  more 
important  than  the  total  noise  power  accompanying  the 
signal. 

A  true  signal  recovery  system  must  be  capable  of 
responding  to  the  signal  'buried  in  noise'.  Moreover,  it 
must  be  able  to  determine  the  amplitude  and  phase 
variations  over  the  desired  frequency  range.  Fig  2  shows 
the  characteristics  of  a  signal  recovery  system  to  be  used 
in  transfer   function  measurements. 
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Phase  Sensitive  Detectors  (PSD)  are  commonly  used  in 
signal  recovery  applications.  A  phase  sensitive  detector 
provides  an  output  proportional  to  the  amplitude  of  the 
signal  and  the  phase  difference  between  the  signal  and  the 
reference.  In  transfer  function  measurements,  a  PSD  based 
system  must  lock  on  to  the  signal  of  interest  over  the 
desired  frequency  range.  Systems  operating  on  the  PSD 
principle   are    termed   lock-in    systems. 

2.1   Principles  of   operation: 

The  requirements  of  a  basic  lock-in  system  are  shown  in 
fig  3.  Lock-in  systems  invariably  use  a  reference  signal 
derived  from  the  excitation  source.  The  correlation 
between  the  signal  of  interest  and  the  reference  signal  is 
tested  by  multiplying  together  the  two  inputs  to  form  the 
product 

Vp(t)    =    r(t)     (    s(t)    +  n(t)    ) 
where 

r(t)   — >   reference  signal 
s(t)   — >   signal  of  interest 
n(t)   — >  noise  and  interference 
The  higher  products  of  multiplication  are  suppressed  by 
a  lowpass  filter  at  the  system  output.  The  final  output 
will  be  a  constant  voltage  proportional  to  the  signal 
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level  when  the  signal  of  interest  and  the  reference  signal 
are  closely  correlated.  The  product  of  n(t)  and  r(t) 
averaged  over  a  period  of  time  will  be  zero. 

2.2  Signal    channel: 

The  signal  channel  mainly  consists  of  the  experimental 
system  driven  by  a  sinusoidal  source.  Commercial  lock-in 
systems  are  usually  provided  with  optional  preamplifiers 
not  merely  to  increase  the  gain  but  also  to  provide  a 
noise  match  to  the  signal  source.  Signal  conditioning 
filters  are  sometimes  provided  to  increase  system 
sensitivity.  There  is  no  need  to  clean  the  signal  prior  to 
detection  in  a  synchronous  system.  The  introduction  of 
filters  in  the  signal  channel  poses  problems  in  phase 
measurements.  The  phase  shifts  introduced  by  these  filters 
have  to  be  compensated  in  the  reference  channel. 
Therefore,    filtering  has   to   be    kept   to   a  minimum. 

2.3  Reference   channel: 

The  main  purpose  of  the  reference  channel  is  to  provide 
a  signal  synchronous  with  the  signal  applied  to  the 
experimental  system.  The  reference  signal  may  be  either  a 
square  wave  or  a  sine  wave  determined  by  the  design  of  the 
phase   sensitive   detector. 


A  phase  shifter  is  provided  in  the  reference  channel  to 
control  the  relative  phase  of  the  two  inputs  to  the 
multiplier.  The  phase  shifter  forms  a  major  component  in 
transfer  function  measurements.  It  has  to  provide  accurate 
phase  shifts  over  the  desired  frequency  range.  The  phase 
shifters  available  are  of  the  voltage  control  type 
requiring  the  tuning  of  either  a  resistor  or  capacitor  to 
obtain  the  proper  phase  shift.  Several  phase  shifter 
circuits  using  transistors,  operational  amplifiers,  Phase 
locked  loops  are  available  in  the  literature.  In  precision 
phase  measurements,  a  calibration  procedure  is  necessary  to 
estimate  the  phase  shift. 

2.4  Phase  sensitive  detector: 

The  multiplier  and  the  lowpass  filter  together 
constitute  the  phase  sensitive  detector.  The  phase 
sensitive  detectors  usually  used  are  of  the  switching 
multiplier  type  as  shown  in  fig  4  due  to  its  exceptionally 
wide  dynamic  range  and  operational  simplicity.  The 
reference  signal  is  a  square  wave  which  controls  an 
electronic  switch.  The  output  of  the  switch  is  passed 
through  a  lowpass  filter  which  provides  a  dc  voltage 
proportional    to  the  amplitude   of   the   signal. 

Fig  5  shows  the  waveforms  at  the  output  of  the  phase 
sensitive    detector    for     different    phase     relations    between 
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Fig   5.  Waveforms  in   a   phase  sensitive   detector 
for  different  phase   relations 

(Meade  M.L,  Lock-in  amplifiers:principles  and  applications,  1983,  p 3 2  —  3 3 
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the  inputs  of  the  multiplier.  The  output  depends  not  only 
on  the  amplitude  of  the  signal  but  also  on  the  phase 
difference  between  the  signal  and  the  reference.  The 
response  of  the  lock-in  system  is  maximum  when  the  inputs 
to  the  multiplier  are  precisely  in  phase.  A  mathematical 
analysis  of  the  lock-in  system  for  both  square  and  sinewave 
references   is   given   in  Appendix  A. 

The  output  lowpass  filter  provides  a  major  means  of 
improving  the  signal  to  noise  ratio  in  lock-in  systems.  The 
purpose  is  to  filter  the  higher  products  of  multiplication 
at  the  final  output.  The  bandwidth  of  the  filter  is 
typically  around  1  Hz.  The  settling  time  of  the  filter 
(inversely  proportional  to  the  bandwidth)  is  important 
otherwise  the  response  will  be  sluggish.  The  lowpass  filter 
may  be  a  single  section  or  double  section  RC  filter  with  a 
roll  off  of  6  dB/Octave  or  12  dB/Octave  respectively  beyond 
the  cut   off   frequency. 

2.5  A  Practical   lock-in   system: 

The  Function  Generator  provides  the  sinusoidal  input  to 
the  lock-in  system.  Most  of  the  function  generators 
available  have  the  'Sync  output1  capability.  In  addition  to 
producing  the  desired  signal,  they  also  produce  a  square 
wave  synchronous  with  the  zero  crossing  of  the  signal.  The 
levels    of    the    square  wave    depend   on   the    termination    on    the 
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'Sync  output'  line.  The  HP3325A  is  one  such  Function 
Generator  and  can  be  used  to  provide  the  input  and 
reference   signals   for   the  lock-in   system. 

The  switching  multiplier  was  implemented  using  the  CMOS 
Analog  Switch  DG305  (SPDT)  as  shown  in  fig  6.  This 
implementation  requires  inverting  and  non-inverting  unity 
gain  amplifiers  in  the  signal  channel  which  causes  errors 
in  precision  phase  measurements.  In  addition,  the  switch 
does  not  work  properly  at  high  frequencies  or  small  signal 
amplitude. 

An  alternative  implementation  is  to  use  a  precision  IC 
Multiplier.  The  AD534  is  a  four  quadrant  analog  multiplier 
without  any  external  trimming  for  offset  and  gain  accuracy. 
The   generalized  transfer   function   is   given   by 

(XI   -  X2)  (Yl   -  Y2) 
SF 


V 


out 


where 


-    (Zl   -   Z2) 


A  =  open  loop  gain  of  output  amplifier,  typically  7  0 
dB  at  dc. 

X,Y,Z   =  differential  input  voltages,  high  impedance 

SF  =  Scale  factor,  pretrimmed  to  10  V. 

The  operation  of  the  AD534  as  a  multiplier  is  described 

by  the  equation 

(XI  -  X2) (Yl  -  Y2)  =  10  V(Z1  -  Z2) 

The    small    signal    bandwidth    is    typically    around    1    MHz    with 

offset  voltage   +   5   mV.    The   lock-in   system   using  an  AD534    is 

shown   in  fig   7. 
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3.    GAIN    PHASE    MEASUREMENTS 

Frequency  response  measurements  refer  to  the 
performance  evaluation  of  a  test  circuit  in  amplitude  and 
phase  when  it  is  subjected  to  a  time  varying  input  voltage. 
The  distortion  in  the  input  signal  has  to  be  minimized  for 
accurate  measurements. The  analog  input  voltage  usually 
provided  in  these  measurements  is  a  sinusoidal  waveform. 
This  is  because  of  availability  of  sinusoidal  waveform 
generators  with  low  distortion  and  the  sinewave  allows  the 
performance   to   be   specified  at    discrete  frequencies. 

3.1  Measurement   procedures: 

Two  measurement   procedures  will   be   discussed. 

1.  Simultaneous  sampling  method. 

2.  Using  lock-in  systems. 

3.1.1  Simultaneous  sampling  method: 

This  method  consists  of  simultaneously  sampling  the 
input  and  output  signals  of  a  test  circuit  for  a  particular 
frequency.  For  proper  measurements,  at  least  two  cycles  of 
the  signals  have  to  be  sampled.  The  gain  and  phase-shift 
of  the  circuit  for  a  particular  frequency  is  obtained  in 
software  based  on  these  samples.  The  maximum  amplitudes  for 
both  the  signals  are  determined  and  the  ratio  of  the 
maximum  output  amplitude  to  the  maximum  input  amplitude 
gives    the    gain    of    the    circuit    for    that    frequency.     The 
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phase-shift  introduced  is  estimated  using  the  difference  in 
zero  crossing  of  the  two  signals  in  a  particular  direction. 
A  knowledge  of  the  sampling  frequency  is  necessary  to 
obtain  the  magnitude  of  the  phase-shift.  The  sign  of  the 
phase-shift  depends  on  which  one  of  these  signals  reaches 
the  positive  peak  first.  However,  it  becomes  difficult  to 
interpret  the  phase  relationship  if  both  the  signals  attain 
the  peak  in  the  same  sampling  interval.  Also,  the  phase 
measurement  is  sensitive  to  any  offset  introduced  by  the 
circuit.  The  accuracy  of  gain  and  phase  measurements  using 
this  method  is  entirely  dependent  on  the  sampling 
frequency.  The  higher  the  sampling  frequency,  more  accurate 
are  the  measurements  .  A  program  to  obtain  simultaneous 
sampling  of  input  and  output  signals  using  the  Keithley 
194A  is   provided   in  Appendix  C. 

3.1.2  Using  lock-in  systems: 

The  key  to  the  measurement  procedure  using  lock-in 
systems  is  to  obtain  the  maximum  average  output  signal.  The 
following  methods  are  possible. 

3.1.2.1  Reference  channel  variations: 

For  any  frequency  setting,  vary  the  phase  of  the 
reference  channel  from  0  to  360  degrees  in  steps  depending 
on    the    accuracy     of    measurements     required.     Obtain    the 
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average  output  voltages  for  each  phase  setting.  From  these 
readings,  the  maximum  average  output  voltage  and  the 
corresponding  phase-shift  can  be  obtained.  If  the  phase- 
shift  obtained  is  greater  than  180  degrees,  the  output 
signal  leads  the  input  signal.  If  the  phase-shift  is  less 
than  or  equal  to  180  degrees,  the  output  lags  the  input. 
The  time  taken  by  this  method  to  perform  a  frequency 
response  measurement  for  a  reasonable  phase  accuracy  of  ±1 
degree  is  quite  large,  (at  least  an  hour!).  The  offset 
introduced  by  the  circuit  affects  the  gain  measurements  but 
not   the  phase  measurements. 

3.1.2.2  Null  shift  procedure: 

The  phase  sensitive  detector  response  is  proportional  to 
Cos0,  where  0  represents  the  relative  phase  of  the  signal 
and  reference  at  the  phase  sensitive  detector  input. 
(Appendix  A) .  The  variation  at  the  output  of  the  phase 
sensitive  detector  as  the  reference  channel  phase  is  varied 
from  0  to  360  degrees  is  shown  in  fig  8.  The  cosine  nature 
of  the  output  can  be  exploited  to  obtain  the  phase-shift 
introduced   by   the    circuit.    The   procedure   is  outlined   below. 

-  Start   from  an  arbitrary    initial    phase   condition    (0   deg) . 

-  The    phase       sensitive    detector    output    is    nulled    by 

adjusting    the    reference    phase,     (quadrature    phase).     This 

establishes   a  quadrature   condition  at   the    phase    sensitive 

detector. 
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-  The  phase-shift   introduced  by  the  circuit  is  then 

established  by  shifting  the  quadrature  phase  by  90 

degrees. 

Frequency  response  measurements  done  using  this  method  are 

susceptible  to  any  offset  introduced  by  the  lock-in  system. 

3.2  Implementation: 

3.2.1  Computerized  Test  system: 

The  basis  of  frequency  response  measurements  using  lock- 
in  systems  is  to  precisely  adjust  the  phase  shifter  in  the 
reference  channel  until  the  inputs  to  the  phase  sensitive 
detector  are  in  phase.  The  phase  shifter  circuits  are 
frequency  dependent  and  physical  tuning  of  circuit  elements 
is  undesirable  in  a  computer  controlled  lock-in  system.  In 
order  to  automate  these  measurements,  function  generators 
with  phase  shifters  and  computer  control  capability  are 
necessary. 

The  block  diagram  of  a  computerized  test  system  is  shown 
in  fig  9.  The  Synthesizer  and  the  Digital  Voltmeter  can  be 
programmed  by  the  controller  HP9  836  over  the  IEEE  4  88  bus. 
A  Printer  and  Plotter  are  also  connected  to  represent  data 
and  results. 

The  HP8904A  Multifunction  Synthesizer  (Option  002)  is  a 
dual  channel  synthesizer  with  the  channels  programmed 
individually.  The  relevant  specifications  of  the 
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synthesizer  are  provided  in  Appendix  B.  Channels  A  &  B  on 
the  synthesizer  are  used  as  the  signal  and  reference 
channels  of  the  lock-in  system  respectively.  The  phase  of 
the  signals  in  both  the  channels  can  be  programmed  from  0 
to  359.9    degrees   in   steps   of    0.1    degree. 

The  Keithley  194A  High  Speed  Voltmeter  is  a  dual  channel 
digital  voltmeter.  The  channels  can  be  programmed 
individually  to  obtain  the  samples  of  the  signal  and  store 
them  in  a  buffer.  Several  mathematical  functions  such  as 
Average,  TRMS,  Peak  etc. ,  are  available  which  are  computed 
based  on  the  samples  obtained.  The  sampling  rate  and  the 
number  of  samples  can  be  programmed  .  More  details 
regarding  the  architecture  and  dual  channel  operation  of 
Keithley  194A  are   provided   in  Appendix   C. 

The  frequency  range  desired  is  upto  100  KHz.  The 
reference  signal  when  using  the  IC  Multiplier  as  the  phase 
sensitive  detector  can  be  either  a  square  or  a  sine  wave. 
The  upper  limit  on  the  frequency  range  is  decided  based  on 
the  reference  signal  chosen.  This  is  due  to  the  maximum 
frequency  specifications  for  sine  and  square  waves  on  the 
HP8904A  synthesizer  (Appendix  B).  Accordingly,  the 
measurements  are  carried  out  from  10  Hz  to  100  KHz  for 
sinewave  reference  and  from  10  Hz  to  50  KHz  for  a  square 
wave   reference. 
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3.2.2  Modifications  to  the  lock-in  system: 

The  sensitivity  of  the  system  is  diminished  considerably 
due  to  the  attenuation  by  10  at  the  multiplier  output.  This 
is  due  to  the  scale  factor,  SFf  being  10.  The  SF  can  be 
trimmed  down  to  3  V.  From  the  manufacturer's 
specifications,  the  maximum  input  allowed  is  1.25  SF. 
Decreasing  the  SF  would  not  help  a  great  deal  in  increasing 
the  sensitivity.  A  much  lower  scaling  voltage  can  be 
achieved  without  any  reduction  of  input  signal  range  using 
a  feedback  attenuator.  The  elimination  of  scale  factor  can 
be  done  in  software  by  multiplying  every  reading  by  the 
scale  factor.  Hardware  scale  factor  elimination  improves 
the  range  in  which  the  voltmeter  has  to  operate.  However, 
two   problems   are  associated  with   this: 

1.  The  bandwidth  of  the  multiplier  is  reduced  to  about  80 
KHz  in  the  presence  of  the  peaking  capacitor.  This  does 
affect   the  performance  when  a   sine  wave   reference   is   used. 

2.  The  output  offset  voltage  increases  by  a  factor  of  10. 
The  offset  affects  the  determination  of  the  quadrature 
phase  especially  when  the  signal  strength  and  the  offset 
are  comparable.  The  error  is  around  10  to  20  degrees.  The 
effect  of  offset  on  phase  determination  is  shown  in  fig  10. 
The  offset  can  be  removed  in  software  as  well  as  hardware. 
The  offset  is  subtracted  from  every  reading  in  software. 
Both   the   techniques   require   the   estimation   of    offset   prior 
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to  the  measurement  process.  Fig  11  shows  the  modified 
multiplier  circuitry  with  unity  SF  and  offset  correction. 
In  our  implementation,  the  SF  and  offset  are  eliminated  in 
hardware  and  software  respectively. 

3.2.3  Algorithm: 

The  null  shift  procedure  outlined  in  section  3.1.2.2 
with  either  a  square  or  sine  wave  reference  signal  is  used 
in  this  implementation. 

An  important  criterion  to  be  considered  is  to  identify 
whether  the  output  signal  leads  or  lags  the  input  signal. 
The  relationship  between  the  reference  signal,  input  and 
output  signals  of  the  test  circuit  for  the  four  general 
phase  shifts  is  shown  in  fig  12.  It  is  apparent  that  when 
the  reference  channel  is  shifted  by  more  than  180  degrees, 
the  output  signal  leads  the  input  or  else  it  lags. 

The  offset  introduced  by  the  lock-in  system  is  removed 
in  software  by  using  the  ZERO  feature  on  the  Keithley  194A 
Voltmeter. 

The  algorithm  for  frequency  response  measurements  using 
the  above  test  system  is  shown  in  the  flowcharts  given  in 
figs  13  &  14.  The  program  written  in  HP  BASIC  3.0,  using 
the  above  algorithm,  is  provided  in  Appendix  D. 
Autoset  routine: 

This  routine  incorporates  both  autophase  selection  and 

autosensitivity.  The  quadrature  phase  occurs  at  two  phases 
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Fig   13.  Main  flowchart  for  transfer  function 
measurements  using   lock-in   systems 
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Fig   1  4.  Autoset  routine 
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180  degrees  apart.  The  quadrature  phase  corresponding  to 
the  positive  zero  crossing  of  the  output  is  determined 
using  the  'Bisection  Search  Method'.  In  situations  wherein 
the  quadrature  phase  cannot  be  determined  properly,  the 
sensitivity  of  the  system  is  increased  by  varying  the 
amplitudes  on   both   the   signal   and   reference    channels. 

3.3   Frequency    response   of   test   circuits: 

Gain  and  phase  measurements  using  the  procedure  outlined 
in  section  3.2.3  were  carried  out  on  the  following  test 
circuits. 

1.  RC  lowpass  filter 

2.  RC  highpass  filter 
3.3.1  RC  lowpass  filter: 

The  transfer  function  of  a  RC  lowpass  filter  is  given  by 

J_o_        m     ___! 

Vi  1  +  jwRC 

With    R  =    22    K    and   C    =    0.022    uF,     the    cut    off    frequency    is 

328.83    Hz.    The   amplitude  and   phase   responses   are    determined 

using   both    sine    wave    and    square   wave    references.    These 

responses     are      compared  with      the   ideal   obtained     from   the 

transfer    function.    Figs    15    &    16    shows    the    amplitude    and 

phase    responses    of    a    RC    lowpass    filter.     The    phase    errors 

with      sine     wave  and   square  wave    references   are   illustrated 

in  fig   17. 
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The  phase  error  tends  to  increase  in  the  high  frequency 
range  when  the  signal  amplitude  is  very  small.  The  error  is 
less  when  a   sine  wave   is   used  as  the   reference. 

3.3.2   RC  highpass   filter: 

The   transfer    function  of    a  RC  highpass   filter    is   given   by 
VQ  jwRC_ 

Vi  1   +    jwRC 

With  R  =  22  K  and  C  =  0.022  uF,  the  cut  off  frequency  is 
328.83  Hz.  The  amplitude  and  phase  responses  of  the 
highpass  filter  are  shown  in  figs  18  &  19.  The  phase  errors 
are    shown   in  fig   20. 

The  phase  error  is  again  large  in  the  high  frequency 
range  where   the   signal   strength   is   large. 

The  amplitude  response  obtained  using  lock-in  systems  is 
a  function  of  the  phase  response  due  to  the  nature  of  the 
measurement  process.  As  evident  from  the  plots,  the  gain 
errors   are  very   small    compared  to   the   phase   errors. 

3.4    Probable   error   sources   in  phase  measurements: 

The  principal  sources  of  error  affecting  phase 
measurements  are  the   following. 

1.  Signal   and  applied   reference   not    strictly    in   phase. 

2.  Errors   due    to  oscillator   distortion. 

3.  Offset  and   drift   in   phase   sensitive    detector. 

4.  Component    sensitivity. 
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3.4.1  In-phase   signal   and   reference: 

The  applied  signal  and  the  derived  reference  have  to  be 
strictly  in  phase  for  accurate  phase  measurements.  Trigger 
errors  may  occur  when  the  reference  is  derived  from  the 
signal.  In  our  case,  the  two  channels  are  independent  and 
according  to  the  manufacturer' s  specifications,  the  phase- 
shift  between  the  two  channels  is  within  ±  0.1  degree  for 
sinewave  signals.  The  phase-shifts  between  the  two  channels 
for  different  signals  were  tested  using  a  dedicated  Gain- 
Phase  meter  HP3  575A.  The  results  of  the  test  are  tabulated 
in  Tables   1   &    2. 

The  phase-shifts  for  equal  amplitude  sinewaves  was  found 
to  be  well  within  the  specifications.  For  unequal 
amplitudes,  the  phase-shift  increases  as  the  frequency 
increases.  The  phase-shifts  for  sine  and  square  waves  were 
found  to  be  large  at  very  low  and  very  high  frequencies. 
The  maximum  phase-shift  was  as  high  as    6   degrees. 

The  phase-shift  between  the  signal  and  the  reference 
channel  on  the  synthesizer  was  added  to  the  phase-shift 
determined  earlier  for  the  highpass  and  lowpass  filters. 
The  resultant  errors  in  the  phase  responses  for  both  the 
filters    are    shown    in    figs    21    &     22. 

3.4.2  Errors   due   to  oscillator   distortion: 

The    response    of    the    phase    sensitive    detector    to    a 

sinusoidal    input    is    cosine    in    nature.     If    the    input    signal 

is    not    strictly    sinusoidal,     the    output    may    not    follow    the 
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Table  1 

Phase  shift  between  channels  on  the  HP8904A 
for  sinusoidal  signals 


1  Frequency 

Phase  sh 

Lft  (Deg)      | 

(Hz) 

1  Ch  A:  1 

% 

1  Ch  A:  5  Vnlc  | 
1  ChB:  1  v£k  | 

1  Ch  B:  1 

1   10 

1      0.0 

1     0.0       | 

1   20 

!     0.0 

1     0.0       | 

1   30 

1      0.0 

1     0.0       | 

1   40 

1      0.0 

1     0.0       | 

1   50 

1      0.0 

1     0.0       | 

I   60 

0.0 

I     0.0       | 

1   70 

0.0 

1     0.0       I 

1   80 

!     0.0 

0.0       | 

1   90 

0.0 

0.0       I 

1   100 

0.0 

0.0       | 

1   200 

0.0 

-0.1       I 

1   300 

0.0 

-0.1       | 

1   400 

0.0 

-0.1       I 

1   500 

0.0 

-0.1       | 

1   600 

0.0 

-0.1       I 

1   700 

0.0 

-0.1      I 

I   800 

0.0 

-0.1       I 

1   900 

0.0 

-0.1       | 

1   1000 

0.0 

-0.1       I 

1   2000 

0.0 

-0.2       I 

1   3000 

0.0 

-0.2       I 

1   4000 

0.0 

-0.2       I 

I   5000 

0.0 

-0.2       I 

1   6000 

0.0 

-0.2       I 

1   7000 

0.0 

-0.3       I 

1   8000 

0.0 

-0.3       I 

1   9000 

0.0 

-0.3       I 

1   10000 

0.0 

-0.3       I 

1   20000     I 

0.0 

-0.4       I 

1   30000     I 

0.0 

-0.5       | 

1   40000     I 

0.0 

-0.6       I 

1   50000     I 

0.1 

-0.7       | 

1   60000     I 

0.1 

-0.7       I 

1   70000     | 

0.1 

-0.8       | 

I   80000     I 

0.1 

-0.8       I 

1   90000     I 

0.1 

-0.8       | 

1   100000    I 

0.2 

-0.9       1 

Note:  Phase  shifts  are  measured  with  respect  to  channel  A 
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Table  2 

Phase  shift  between  channels  on  the  HP8904A 
for  sine  and  square  waves 


1  Frequency 

Phase  shift  (Degrees) 

1 

Cr 

i  A:  Sine 

wave 

I   Ch  A:  Square 

wave    I 

(Hz) 

a 

i  B  :  Square 

wave 

I   Ch  B:  Sine 

wave      I 

A: 

}   VP*  ! 

A 

:  ?  Vpk 

1   A: 

j  vpk  ! 

A 

:  I   VPk  ! 

B: 

1  vpk  1 

B 

5  1  vpk 

1   B: 

1  vpk  i 

B 

:  5  vpk  1 

1   10 

-0.8    | 

-0.9 

-0.9    I 

-0.9    1 

1   20 

-0.4    I 

-0.4 

-0.4    | 

-0.4    | 

1   30 

-0.2    I 

-0.3 

-0.3    1 

-0.3    I 

I   40 

-0.1    1 

-0.2 

-0.1    I 

-0.2    | 

1   50 

-0.1    I 

-0.1 

-0.1    1 

-0.1    I 

1   60 

0.0    I 

-0.1 

-0.1    I 

-0.1    | 

I   70 

0.0    I 

-0.1 

0.0    I 

-0.1    I 

1   80 

0.0    I 

-0.1 

0.0    | 

-0.1    I 

1   90 

0.0    | 

0.0 

0.0    I 

0.0    I 

1   100 

0.0    | 

0.0 

0.0    | 

0.0    | 

1   200 

0.1    I 

0.0 

0.1    | 

0.0    I 

1   300 

0.1    I 

0.0 

0.1    I 

0.0    | 

1   400 

0.1    I 

0.0 

0.1    | 

0.0    I 

1   500 

0.1    | 

0.0 

0.1    I 

0.0    | 

1   600 

0.1    | 

0.0 

0.1    I 

0.0    | 

1   700 

0.1    I 

0.0 

0.1    I 

0.0    | 

1   800 

0.1    | 

0.0 

0.1    | 

0.0    1 

1   900 

0.1    I 

0.0 

0.1    I 

0.0    I 

I   1000 

0.1    1 

0.0 

0.1    I 

0.0    I 

I   2000 

0.1    1 

-0.1 

0.0    1 

-0.1    I 

I   3000 

-0.1    1 

-0.3 

-0.1    1 

-0.3    1 

1   4000 

-0.2    | 

-0.5 

-0.3    1 

-0.5    | 

1   5000 

-0.4    I 

-0.6 

-0.4    I 

-0.6    I 

I   6000 

-0.5    I 

-0.8 

-0.6    I 

-0.8    | 

1   7000 

-0.7    I 

-0.9 

-0.7    1 

-0.9    I 

1   8000 

-0.8    I 

-1.1 

-0.8    I 

-1.2    I 

1   9000 

-0.9    I 

-1.2 

-1.0    I 

-1.3    I 

1   10000 

-1.1    1 

-1.4 

-1.2    I 

-1.5    I 

I   20000 

-2.5    I 

-2.9 

-2.6    I 

-3.0    I 

I   30000 

-3.6    | 

-4.0 

-3.8    | 

-4.3    | 

1   40000 

-4.5    I 

-5.0 

-4.9    I 

-5.5    I 

1   50000 

-5.6    I 

-6.0 

-5.9    1 

-6.6    | 

Note:  Phase  shifts  are  measured  with  respect  to  sine  wave 
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Cos8  law.  It  is  necessary  for  precision  phase  measurements, 
the  distortion  of  the  input  signal  be  very  small. 

Oscillator  distortion  results  in  harmonics.  The  effect 
of  harmonic  components  is  to  shift  the  zero  crossings  of 
the  signal  relative  to  its  fundamental  component.  As  a 
result,  there  exists  an  arbitrary  phase  relationship 
between  the  input  and  the  reference  signal  causing  errors 
in  measurements. 

3.4.3  Offset  and  drift  in  phase  sensitive  detector: 

The  offset  voltage  is  a  main  source  of  error  in  phase 
measurements  based  on  the  zero  crossing  of  the  signals.  The 
effect  of  offset  was  shown  in  section  3.2. 

The  output  drift  in  the  phase  sensitive  detector  and 
phase  drift  in  the  reference  channel  can  cause  errors  when 
the  reference  phase  is  being  adjusted  for  a  null  output. 

3.4.4  Component   sensitivity: 

The   measurements    carried    out    have    to    be    compared    if 

possible  with    the    ideal    results    in    order    to    determine    the 

phase   error.    The   ideal    results   are   calculated   based   on   the 

transfer    function    of    the    test     circuit.     The     circuit 

components    that    are    actually    used    have    to    be    measured 

precisely    in    order    to    determine    the    errors.     The    resistor 

and   capacitor    for    the   highpass   and   lowpass    filters   were 

supposed   to    be   22   K  and   0.022   uF.    The  measured  values   were 

22. 3K    and    0.017    uF    respectively.    This    corresponds    to   a    22% 
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decrease  in  the  RC  product.  The  phase  errors  for  different 
percentage  changes  in  RC  product  for  both  highpass  and 
lowpass  filters  is  shown  in  figs  23  &  24.  It  is  apparent 
that  the  phase  errors  are  dependent  on  frequency,  circuit 
components  and  function  of   the   circuit. 
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4.  SOFTWARE  SIMULATION  OF  GAIN  PHASE  MEASUREMENTS 

The  determination  of  amplitude  and  phase  responses  of 
linear  systems  is  based  on  the  simultaneous  sampling  of 
input  and  output  signals  of  the  system.  Since  the  input  and 
output  signals  are  sinusoidal,  sine  wave  curve  fitting 
routines   can   be   used  to  obtain  the   gain  and   phase-shift. 

The  sine  wave  curve  fitting  routines  provide  estimates 
of  the  amplitude,  frequency,  offset  and  phase  to  a  digital 
data  record.  Several  curve  fitting  routines,  both  closed 
form  and  iterative,  are  available  in  the  literature.  They 
include  the  closed  form  algorithm  [2] ,  the  three-parameter 
(known  frequency)  least  squares  fit  algorithm  and  the  four- 
parameter  least  squares  fit  algorithm  [3].  In  our 
application,  since  the  input  frequency  is  known,  the  three 
parameter   least   squares  fit  algorithm  would  be   ideal. 

An  alternative  approach  would  be  to  test  the  correlation 
between  the  input  and  output  signals  of  the  system.  The 
basis  of  transfer  function  measurements  using  lock-in 
systems  was  to  test  the  correlation  of  the  signals  at  the 
input   of   the  multiplier. 

A  brief  dicussion  of  the  above  algorithms  is  provided  in 
the  next  section  followed  by  the  implementation  of  the 
algorithms  on  two  test   circuits. 
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4.1   Three  parameter   least   squares   fit  algorithm: 

The  three  parameter  algorithm  estimates  the  amplitude, 
offset  and  phase  of  the  digital  data  record  for  known 
frequency.  The  algorithm  provides  a  closed  form  solution 
and  the  accuracy  of  the  estimations  will  be  poor  if  the 
actual  frequency  of  the  signal  sampled  differs  from  the 
frequency  used  in  the  algorithm. 
Algorithm: 

A  detailed  description  of  the  algorithm  along  with  its 
derivation  may  be  found  in  [3] .  The  input  to  the  algorithm 
consists  of  a  data  record  of  N  samples  yn  taken  at  times 
tn,   where   n  =   0,1,... N-l. 

The  algorithm  provides   the   solution   in  the   form 
yn'    =   A  Cos(wtn)    +   B   Sin(wtn)    +    C 
where 

w     — >      known  angular    input   frequency 
t     — >      sample   times 
However,    we    are    interested    in    obtaining    a    solution    of    the 
form 

yn'    =   R  Sin(wtn  +    0)    +    C 
The    amplitude    R    and    the    phase    0    of    the    sinewave    can    be 
obtained   by    using   the   following   relations 

R  =    (A2  +    B2)1//2 
0    =    tan   _1 (A/B) 
The  offset   present    in   the   data    record   is   estimated   by   C. 
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The  amplitude  and  offset  are  estimated  properly.  The  input 
signal  phase  range  is  from  0  to  360  degrees,  while  the 
estimated  phase  does  not  exceed  +  90  degrees.  A  suitable 
phase  correction  scheme  has  to  be  applied  to  the  algorithm 
to  obtain  the  proper  phase  shift.  In  applications  such  as 
transfer  function  measurements,  where  we  need  to  estimate 
the  phase-shift  between  two  signals,  it  is  necessary  that 
the  phase  be  defined  from  a  common  reference.  It  is 
convenient  to  determine  the  phase  using  the  positive  B  axis 
as   the   reference. 

The  phase  correction  scheme  depends  on  the  signs  of  A 
and  B.    The  actual    phase  0   may   be  obtained  as   follows: 

0      =      9  — >   1st  quadrant 

0      =    it    +    0        — >   2nd  quadrant 

0      =     TT   +    6        — >   3rd  quadrant 

0      =    2.TJ  +    0        — >    4th  quadrant 
The  above   scheme   is   illustrated   in  fig   25. 

The  gain  and  phase  measurements  for  a  particular 
frequency  using  the  three  parameter  least  squares  fit 
method   is    implemented  as   follows: 

1.  Inputs:    input   samples [yin] ,    output   samples [yout] , 

number    of    sampl es [num] ,    sample    times[xn]    and 
f  requency[f  req]  . 

2.  A    least    squares    fit    of    the    input     samples    provides 
estimates   of   the   input   amplitude,    offset  and  phase. 
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Ill  Quadrant 


IV  Quadrant 


Fig   25.  Phase   correction   scheme   used   in   three 
parameter  least  squares  fit  algorithm 
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3.    A    least    squares    fit    of    the    output     samples    provides 
estimates   of   the  output   amplitude,    offset  and   phase. 

output   amplitude 


4.  Gain(dB)    =   20    log 

input   amplitude 

5.  Phase-shift  =  output   phase  -   input   phase 


4.2  Cross  correlation  using  Past  Fourier  Transforms: 

Correlation  tests  the  similarity  of  two  data  sets.  The 
correlation  between  input  and  output  samples  is  tested  by 
comparing  them  both  directly  superposed.  The  correlation 
will  be  large  at  some  value  of  time  t  if  the  output  signal 
lags  the  input  signal.  Likewise,  the  correlation  will  be 
large  for  some  negative  value  of  t  if  the  output  leads  the 
input. 

The  correlation  between  the  input  and  output  samples  is 
computed  as  follows: 

1.  Obtain  the  Discrete  Fourier  Transform  (DFT)  of  the  two 
data  sets. 

2.  Multiply  one  resulting  transform  by  the  complex 
conjugate  of  the  other. 

3.  Obtain  the  Inverse  Discrete  Fourier  Transform  (IDFT)  of 
the  product. 

The  result  (vector  rk)  will  formally  be  a  complex  vector  of 
the  sample  length  N.  The  correlation  at  zero  lag  is  in  rQ, 
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the  first  component;  the  correlation  at  lag  1  is  in  r*f  the 
second  component;  the  correlation  at  lag  -1  is  in  rN_lf  the 
last  component  etc.  The  gain  can  be  calculated  by  obtaining 
the  maximum  amplitude  of  input  and  output  samples  for  a 
particular  frequency  and  using  the  relation 

max.  output  amplitude 


Gain(dB)  =  20  log   

max.  input  amplitude 


Discrete  Fourier  Transform: [4] [5] 

The  DFT  of  a  function  h^  from  a  finite  number  of  sampled 
points  N  is  given  by 

N-l 

Hn  =7"     hk   e2TT  lkn/N  n   =   -N/2    ••    N/2 

k=0 

The   IDFT   is    given   by 

N-l 


=  _LY~~    Hn   e-2TTil 
nZ — 


hk  -  "■   ikn/N 


N 
n=0 

The  only  differences  in  the  above  two  equations  are 

(i)    Changing  the  sign  in  the  exponential 

(ii)   Dividing  the  result  by  N 

Therefore,  a  single  routine  for  calculating  DFT  can  also, 

with  slight  modifications,  calculate  the  IDFT. 

The  computation  of  DFT  using  the  above  equation  results 

in  redundant  calculations  [4].  This  redundancy  can  be 
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overcome  by  computing  the  DFT  using  an  algorithm  called  the 
Fast  Fourier  Transforms  (FFT) . 

The  computation  of  the  DFT  using  the  FFT  is  based  on  the 
Danielson  Lanczos  Lemma  which  states  that  a  DFT  of  length  N 
can  be  obtained  from  the  two  discrete  fourier  transforms, 
each  of  length  N/2.  One  of  the  two  is  formed  from  the  even 
numbered  points  of  the  original  N,  the  other  from  the  odd 
numbered  points.  The  Danielson  Lanczos  Lemma  can  be  used 
recursively  until  we  have  subdivided  the  data  all  the  way 
down  to  transforms  of  length  1.  The  fourier  transform  of 
length  one  is  just  the  identity  operation  that  copies  its 
one  input  number  into  its  output  slot.  The  output  slot  for 
a  particular  decomposition  is  the  location  obtained  by  bit 
reversal  of  the  original  sample  location.  The  DFT  can  be 
computed  quickly  but  the  number  of  samples  has  to  be  an 
integer  power  of  2. 

The  above  algorithms  have  been  implemented  in  a  ' C 
VAX/VMS  environment  and  the  source  code  listings  are 
provided  in  Appendix  D. 

4.3  Frequency  response  of  test  circuits: 

The  RC  lowpass  and  highpass  filters  (R  =  22. 3K,  C  = 
0.017uF)  used  in  section  3.3  were  used  again  as  the  test 
circuits.  The  input  to  the  test  circuits  was  a  sinewave  of 
amplitude  5  V  k.  The  input  and  output  samples  over  the 
frequency  range  10  Hz  to  100  KHz  for  both  the  circuits  were 
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obtained  using  the  program  provided  in  Appendix  C.  Since 
the  number  of  samples  for  correlation  has  to  be  an  integer 
power  of  2,  the  samples  were  obtained  accordingly.  The 
maximum  sampling  frequency  on  the  Keithley  194A  placed 
restrictions  on  the  number  of  samples  per  cycle.  Table  3 
gives  the  number  of  points  sampled  per  cycle  for  different 
frequencies. 

Table  3 
Number  of  points  sampled  for  different  frequencies 
Frequency  (Hz)  Samples  per  cycle 

10  512 


IK  512 

2K  256 

3K  256 

4K  128 

5K  12  8 

6K  128 

7K  12  8 

8K  64 

9K  64 

10K  64 

2  OK  3  2 

3  OK  3  2 
40K  16 
50K  16 

6  OK  16 

7  OK  8 
80K  8 
90K  8 
100K  8 


The  phase  accuracy  when  using  the  cross  correlation 
method  is  dependent  on  the  number  of  points  sampled  per 
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cycle.  Table  4  shows  the  phase  resolution  (in  degrees)  as  a 

function  of  the  number  of  points  per  cycle. 

Table  4 

Phase  resolution  as  a  function  of  points  sampled  per  cycle 
in  cross  correlation  method 

Points  per  cycle  Phase  resolution 

(degrees) 

1024  +    0.3516 

512  +   0.7031 

256  +    1.4063 

128  ±    2.8125 

64  ±    5.625 

32  ±   11.25 

16  ±    22.5 

8  +    45.0 

The  amplitude  and  phase  responses  of  the  RC  lowpass 
filter  obtained  using  both  the  software  techniques 
implemented  are  shown  in  figs  26  &  27  respectively.  The 
amplitude  response  of  the  filter  using  least  squares  fit 
method  provides  better  results  than  the  correlation  method. 
This  is  because  the  least  squares  fit  provides  an  amplitude 
estimate  for  the  data  record  whereas  the  cross  correlation 
method  uses   the  maximum  amplitude   in  the   data   record. 

The  phase   error    in  the   least   squares  method   is  generally 

less   than   the    cross    correlation  method.    The    phase    error    in 

correlation    increases    when    the    sample    size    decreases     (at 

high    frequencies).    However,     the    phase    error    at    10    KHz    is 

extremely    large    (around    22    degrees) .    This    error    is    of    the 

same  magnitude    using   both    the  methods    implying    sampling 

error   on  part   of    the  Keithley   194A  Voltmeter. 
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The  amplitude  and  phase  responses  of  the  RC  highpass 
filter  are  shown  in  figs  28  &  29  respectively.  The 
observations  are  similar  to  that  made  for  the  lowpass 
filter.  Extremely  large  phase  errors  are  observed  at  30  & 
400  Hz  in  both  the  methods  due  to  sampling  error. 

The  error  in  phase  measurements  for  lowpass  and  highpass 
filters  are  shown  in  figs  30  &  31  respectively. 

The  sampling  resolution  of  the  Keithley  194A  Voltmeter 
and  the  constraint  on  number  of  samples  per  cycle  results 
in  errors  since  a  complete  cycle  of  the  signal  is  not 
sampled. 

The  correlation  method  was  tested  for  the  two  circuits 
with  simulated  samples  of  512  samples  per  cycle  over  the 
frequency  range  1  KHz  to  100  KHz.  The  error  in  phase 
measurements  are  tabulated  in  Table  5.  As  expected  the 
errors  are  well  within  the  maximum  error  (resolution/2) . 

The  least  squares  fit  method  was  tested  for  the  two 
circuits  with  simulated  samples  of  8  samples  per  cycle  over 
the  frequency  range  10  Hz  to  100  KHz.  The  phase  error  for 
both  the  circuits  were  less  than  0.1  degree. 
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Table  5 

Phase  error   obtained  for   lowpass  and  highpass  filters  with 
simulated   samples   using  cross   correlation  method 


Frequency  Phase  error    (degrees) 

(KHz)  RC  lowpass  filter  RC  highpass   filter 

0.27379 

-0.098199 

0.231922 

-0.336538 

-0.122314 

-0.216258 

-0.083437 

0.191506 

-0.141759 

0.294623 

-0.203722 

0.0  986  26 

-0.101795 

-0.222054 

-0.30223 

0.343626 

0.300673 

0.267266 

0.24054 

Phase  error  =   Ideal   value  -   Measured  value 
Phase  resolution  for   512    samples   per   cycle   is     +0.7031   deg. 


1 

-0.27379 

2 

0.098199 

3 

-0.231922 

4 

0.336538 

5 

0.122314 

6 

0.216258 

7 

0.083437 

8 

-0.191506 

9 

0.141759 

10 

-0.294623 

20 

0.203722 

30 

-0.0986  26 

40 

0.101795 

50 

0.222054 

60 

0.30223 

70 

-0.343626 

80 

-0.300673 

90 

-0.267266 

100 

-0.24054 
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5.    CONCLUSIONS 

The  purpose  of  this  work  was  to  compare  the  performance 
of  transfer  function  measurement  schemes  in  hardware  and 
software. 

The  hardware  approach  using  computer  controlled  lock-in 
systems  provides  phase  accuracy  dependent  on  the  resolution 
of  the  phase  shifters  available.  Using  the  precision  IC 
multiplier  AD534  in  the  phase  sensitive  detector  enables  us 
to  use  either  a  square  or  sine  wave  reference.  The  response 
is  generally  better  to  a  sine  wave  reference.  It  was  found 
that  the  phase  errors  were  primarily  due  to  the  phase  shift 
between  the  signal  and  reference  channels  on  the 
synthesizer.  The  method  provides  scope  for  varying  signal 
and  reference  amplitudes,  thereby  increasing  the 
sensitivity  of  the  system.  However,  the  user  must  be  aware 
of  the  maximum  output  amplitude  of  the  test  circuit  as  this 
may  exceed  the  range  of  input  amplitudes  on  the  multiplier. 
The  amplitude  response  is  a  function  of  the  phase  response 
due   to  the  nature  of   the  measurement   process. 

The  software  techniques  require  simultaneous  sampling  of 
the  input  and  output  signals  of  the  test  circuit.  The 
number  of  points  sampled  per  cycle  plays  a  critical  role  in 
phase  estimation  using  the  cross  correlation  method.  There 
should  be  at  least  180  samples  per  cycle  for  a  phase 
accuracy    of  +    1    degree.    The   cross    correlation   method    using 
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Fast  Fourier  Transforms  requires  that  the  number  of  sample 
points  be  an  integer  power  of  2  in  computing  the  Discrete 
Fourier  Transform.  The  computation  of  the  DFT  using 
straight  forward  method  does  not  place  any  restriction  on 
the  sample  size,  but  is  slower.  However,  the  phase 
resolution   still   depends  on   the  number   of   points   per   cycle. 

The  least  squares  fit  method  is  more  robust  and  it 
provides  phase  accuracy  to  within  +  0.1  degrees  with  8 
samples  per  cycle.  The  only  constraint  on  this  method  is 
that  the  frequency  used  in  the  algorithm  be  the  same  as 
that  of   the   signal    sampled. 

The  software  techniques  rely  on  the  accurate 
simultaneous  sampling  of  input  and  output  signals.  The 
sampling  resolution  of  the  Keithley  194A  Voltmeter  and  the 
constraint  on  the  number  of  points  per  cycle  causes 
difficulties  in  sampling  a  complete  cycle.  An  improvement 
to  the  sampling  program  would  be  to  sample  two  cycles  of 
the  signal  and  then  select  samples  of  one  cycle  in 
software.  Care  should  be  taken  to  satisfy  the  constraint  on 
number  of  points  per  cycle.  Controlled  triggering  on  the 
Keithley  194A  may   alleviate   the   problem. 
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APPENDIX  A 
Mathematical  Analysis  of  lock-in  systems 
The  phase  sensitive  detector  of  the  lock-in  system  can 
be  modelled  as  an  ideal  multiplier  followed  by  a  lowpass 
filter.  The  excitation  of  the  experiment  is  sinusoidal  at 
frequency  ws  and  provides  a  reference  which  is  to  be  used 
for  detection  of  the  signal  at  the  output  of  the 
experimental  system.  This  output  signal  is  at  the  same 
frequency  as  the  excitation,  but  suffers  a  phase-shift  0S 
in  the  experiment.  The  reference  is  applied  to  the 
multiplier  via  a  phase  shifter.  The  mathematical 
relationship  for  sine  and  square  wave  references  are 
derived  below.  The  ideal  multiplier  model  for  a  phase 
sensitive   detector   is   shown   in  fig  Al. 

1.  Square  wave  reference: 

The  signal  and  the  reference  voltages  can  be  expressed 

as 

s(t)  =|2  Vs  Cos(wst  +  es) 

r(t)  =_4.Vr  [Cos(wrt+  9r)  -  J^  Cos  3(wrt+  8r) 
TT  3 

+  J.  Cos  5(wrt+6r)  -+..] 
5 


where 


Vc      — >      RMS   value   of    the   signal 

Vr  — >   Peak  value  of  the  reference 


The  product  Vp(t)  is  given  by 
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vp(t)    =  s(t)    *   r(t) 

=  2J2  vs  vr  [  cos(wst  +  wrt  +  es  ±  er) 

-_1   Cos(wst  +    3wrt  +    9S  ±   39r) 
+  _1_  Cos(wst  +    5wrt  +    6S  ±    56r) 
-   +    ..] 
The   signal   and  the   reference  are  at   the   same  frequency.    The 
lowpass   filter    cuts   off  well    below    the    reference    frequency 
which   eliminates   the   higher    products   of  multiplication.    The 
d.c   component   of   the   final    output    is   given  by 

vavg   =   2J2  Vs  Vr   AL(0)    Cos(0s   -   8r) 

TT 

where  AL(0)    is  the   gain   of   the  lowpass  filter   at   dc. 
When    the    signal    and    the    reference    at    the    input    of    the 
multiplier    are    precisely    in    phase     (8S    =    9r)  ,     the    average 
output   voltage   is    given    by    (assuming  AL(0)    and    the   peak 
reference  voltage   to   be   unity) 

vavg  "   2j  vs 

Therefore,  the  scale-factor  to  be  used  to  obtain  the  RMS 

value  of  the  output  signal  of  the  experimental  system  is 
TT 
2/2 
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2.  Sinewave  reference: 

The  signal  and  the  reference  voltages  can  be  expressed  as 


s(t)  =J2  Vs  Cos(wst  +  es) 


r  (t)    =$2  Vr    Cos(wrt  +    9r) 
where  Vs  and  Vr   are   the  RMS  values   of    the    signals. 
The    product  V_(t)    is    given   by 

Vp(t)    =  Vs  Vr    Cos((ws  +    wr)t  +    6S  ±    8r) 
Since    ws    =     wr     and    neglecting    higher     products    of 
multiplication,    the  final    output   is   given  by 

vavg<fc>    "  vs  vr  al(°>    Cos(0s   "   er> 
when  0S   =   9r, 

Vavg(t)    =  Vs  Vr  (assuming  AL(0)    to  be    unity) 

The   signals   on  the  HF8  904A  synthesizer   are  provided   in   peak 

or   peak  to   peak  values.    Assuming  Vr/pkN    to   be   unity, 

VS       =J2       VaVg 

Therefore,     the    scale-factor    to    be    used    to    obtain    the    RMS 
value   of    the   output    signal   of   the  experimental    system   is/2. 
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APPENDIX  B 


HP8904A  Multifunction  Synthesizer  Specifications  [9] 


Wave 
Form 

Frequency 
Range 

Frequency 
Resolution 

Phase 
Range 

Phase 
Resolution 

Amplitude 
Range 

^ 

0  Hz  to  600  kHz 

0.1  Hz 

359.9°    or 
6.282  rad 

0.1°    or 
0.001  rad 

Oto  10  V(p-p) 
(500  load) 

*V 

0  Hz  to  50  kHz 

0.1  Hz 

359.9°    or 
6.282  rad 

0.1°   or      2 
0.001  rad 

Oto  10  V(p-p) 
(50H  load) 

^ 

0  Hz  to  50  kHz 

0.1  Hz 

359.9°   or 
6.282  rad 

0.1°   or      2 
0.001  rad 

Oto  10  V(p-p) 
(500  load) 

Tj 

0  Hz  to  50  kHz 

0.1  Hz 

359.9°   or 
6.282  rad 

0.1°   or      2 
0.001  rad 

Oto  10  V(p-p) 
(50fi  load) 

Gaussian 
Noise 

(N/A) 

(N/A) 

(N/A) 

(N/A) 

0  to  10  V  (p-p) 
(50Q  load) 

(N/A) 

(N/A) 

(N/A) 

(N/A) 

±  10  V 

(open  circuit) 

Table  Bl.  HP8904A  Multifunction  Synthesizer  Capabilities 


Operating  Considerations:  Option  002 

Destination  Control 

Option  002  does  not  allow  you  to  change  the  Destination 
of  the  two  channels.  Channel  A  is  always  designated  for 
Output  1;  Channel  B  is  always  designated  for  Output  2. 

Abbreviated  Specifications:  Option  002 

Output  1  to  Output  2  phase  accuracy  (sine  waves  of  the 
same  frequency) : 

+0.1  degree  or  30  nsf  0.1  Hz  to  100  KHz ,  whichever  is 
greater. 

Bl 


APPENDIX  C 
Dual  Channel  operation  of  Keithley  194A  [8] 

The  Keithley  194A  is  a  high  speed  dual  channel  digital 
voltmeter.  The  channels  can  be  programmed  to  obtain  samples 
of  the  signal  and  store  them  in  a  buffer.  Several 
mathematical  functions  such  as  Average,  TRMSf  Peak  etc.  , 
are  available.  These  functions  are  calculated  based  on  the 
samples  obtained.  The  sampling  rate  and  the  number  of 
samples  can  be  programmed. 

The  architecture  of  the  Keithley  194A  Dual  Voltmeter  is 
shown  in  fig  CI.  Each  channel  has  its  own  A/D  converter, 
which  digitizes  the  analog  input  signal,  and  a  6 4K  byte 
measurement  buffer  for  storage  of  samples.  The  mathematical 
functions  are  calculated  based  on  these  samples  and  placed 
in  a  reading  buffer.  Each  channel  can  process  raw  data 
independently.  However,  the  two  channels  share  the  same 
microcomputer,  reading  buffer  and  the  IEEE  488  bus.  The 
voltmeter  output  always  comes  from  the  reading  buffer 
whether  or  not  that  buffer  has  been  enabled. 
1.  Reading  buffer  disabled:  In  this  mode,  the  reading 
buffer  has  an  effective  size  of  one,  with  both  channels 
sharing  a  same  location.  The  location  will  be  continously 
updated  with  data  from  either  channel  as  the  measurement 
from  that  channel  is  completed,  and  the  reading  is 
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processed.  Any  previous  reading  will  be  overwritten  by  the 
new  one. 

2.  Reading  buffer  enabled:  In  this  case,  a  maximum  of  100 
readings  from  both  channels  can  be  stored  in  the  order  the 
readings  become  available. 

Therefore,  care  has  to  be  taken  to  make  certain  that 
data  sent  over  the  bus  comes  from  the  desired  channel.  The 
following  solutions  are  possible: 

1.  Turn  off  the  unused  channel. 

2.  Discretionary  use  of  trigger  modes  for  both  channels. 

3.  Enable  the  reading  buffer.  Using  suitable  data  formats, 
seperate  the  data  using  the  channel  number  in  the  suffix  in 
software. 

Triggering: 

There  are  two  steps  involved  for  initiating  a 
measurement  sequence.  First,  the  A/D  converter  must  be 
armed,  so  that  it  is  running  and  processing  data.  Secondly, 
the  unit  must  be  triggered  by  the  appropriate  trigger 
stimulus  (determined  by  the  selected  source)  before  it 
performs   a  measurement   sequence. 

Samples  and  Rate  selection  considerations: 

The  sampling  rate  has  to  be  sufficiently  high  so  as  not 
to  lose  important  information  present  in  the  original 
signal.  The  sampling  frequency  must  be  at  least  twice  the 

C3 


highest  frequency  component  in  the  measured  signal. 
Otherwise,  aliasing  results.  The  sampling  duration  or  the 
number  of  samples  is  based  on  how  many  cycles  or  how  much 
of  a  single  cycle  we  wish  to  sample.  The  maximum  number  of 
samples  is  32,768.  The  overall  resolution  and  accuracy  of 
the  measurement  has  to  be  considered  when  selecting  the 
sampling  rate.  The  A/D  converter  operates  with  16  bit 
resolution  for  sampling  rates  less  than  or  equal  to  100 
KHz.  Above  100  KHz,  it  operates  with  8  bit  resolution. 

The  resolution  of  the  programmed  sampling  interval  is 
0.1  usee.  If  the  programmed  sampling  rate  results  in  an 
interval  below  the  resolution,  the  voltmeter  automatically 
adjusts  the  frequency  to  the  nearest  whole  interval.  For 
example,  a  sampling  rate  of  35  KHz  results  in  an  interval 
of  28.571  usee.  The  voltmeter  adjusts  the  interval  to  28.5 
usee  giving  a  sampling  rate  of  35.0877193  KHz. 

Simultaneous  sampling  of  channels: 

The  architecture  of  the  Keithley  194A  poses  problems  for 
simultaneous  sampling  of  channels.  The  channels  are 
programmed  in  the  waveform  mode  (FOX) .  The  channels  can  be 
triggered  simultaneously  by  using  the  following  trigger 
commands. 

Channel  1:  T3X  (Single,  Group  Execute  Trigger) 
Channel  2:  T25X  (Single,  Other  channel) 
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Channel  1  can  also  be  programmed  to  trigger  at  a 
predetermined  slope  and  trigger  level. 

The  samples  are  stored  in  the  respective  measurement  buffer 
in  the  locations  pointed  by  the  Bl  and  B2  pointers.  Once 
the  samples  are  transferred  to  the  controller  over  the  IEEE 
488  bus,  the  channels  have  to  be  programmed  to  some  other 
function  (say  FIX)  in  order  to  cancel  the  waveform  output 
to  enable  further  access  to  samples  from  location  0.  The 
sampling  rate  for  a  particular  frequency  can  be  obtained 
from  the  voltmeter  using  the  status  command  U6.  The  setup 
for  simultaneous  sampling  is  shown  in  fig  C2.  The  program 
for  simultaneous  sampling  is  given  in  page  C7. 
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10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

13  0 

140 

150 

160 

170 

180 

190 

200 

210 

220 


23  0 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
43  0 
440 
450 
500 
510 
520 
53  0 
540 
550 
560 
57  0 


PROGRAM   TO    OBTAIN    SIMULTANEOUS    SAMPLING    OF    CHANNELS 
ON    THE    KEITHLEY    194A 

T.R. RAMESH  FILE:    TFM_DATA 


OPTION   BASE    0 

INTEGER    I, J, Z, Length, Points 

DIM  A(2500) ,B(2500)  fFreq(100)  rSamplesf 


Fn_generator   =  717 
Voltmeter   =7  26 
REMOTE    Fn_ generator 
REMOTE  Voltmeter 
CLEAR  Fn_generator 
CLEAR  Voltmeter 
PRINTER    IS    CRT 


HP- IB 
HP- IB 


Reset 
Reset 


address 
address 


D$[8] ,    S$[8] 

of   HP8904A 

of   Keithley   194A 


HP8904A 
Keithley  194A 


Initialize  HP8904A 


DISP   "Initializing   Fn_ generator. .. ■ 

OUTPUT   Fn_ generator ;"PS"  !    Preset 

OUTPUT   Fn_ generator ;"FC10FW    !    Float    control    Ch.A  off 

OUTPUT  Fn_ generator ;"GM0"         !    Channel    configuration 

!   mode 

OUTPUT    Fn_ generator ;HFRA10 HZ    APA5VL    PHA0DG  WFASI" 

•    Ch.A:    Frequency   =   10    Hz ,    Phase   =    0   Deg 
!    Amplitude  =   5   V   u,    Waveform  =   Sine 

OUTPUT   Fn_ generator ; "0020F    >"    !    Disable   output  B 

Generate   a  frequency   array   from  10    Hz   to  100   KHz    in 
logarthmic  steps 

DATA  10.0,20.0,3  0.0,40.0,50.0,60.0,7  0.0,80.0,90.0 

J   =   0 

FOR  Z    =    0    TO   4 

FOR    I    =   0    TO    8 
READ    Freq(I  +  J) 
Freq(I+J)    =   Freq (I+J) *10"Z 
IF  Freq(I+J)    >100000.0    THEN  GOTO   Start 
Length  =  Length  +  1 
NEXT    I 
J  =   J+9 
RESTORE 
NEXT  Z 


Create  a   data  file   for   storage   of    samples.    The  file 
will   be    stored  on  a   disk   in  the   right  drive   of   the 
HP9836 


CREATE    ASCII    "HPFILTER : INTERNAL, 4" ,100  0 
ASSIGN    @Ofile    TO    "HPFILTER: INTERNAL, 4" 
OUTPUT    @Ofile; Length 
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Points  per   cycle 
Sampling  frequency 


6  00  ! ! 

610  !  Obtain  the  input  and  output  samples  for  all  the    ! 

620  !  frequencies                                     ! 

63  0  ! ! 

700  DISP   "Initializing  Voltmeter. ... " 

710  FOR   I   =  0    TO  Length-1 

720  Points  =   512 

730  Samples  =    Points-1    *   Freq(I) 

7  40  IF   Samples    >    l.E+6    THEN 
750  Points  =   Points/2 
770  GOTO   730 
7  80  END   IF 

790  REDIM  A(Points-l) ,B(Points-l) 

800  OUTPUT    Fn_generator;"FRA"&VAL$(Freq(I) )&"HZn 

810  OUTPUT  Voltmeter ;"C1G1M2R0X"      !    format, autorange 

820  OUTPUT  Voltmeter ; wN0f "&VAL$ (Points) &"S1, "&VAL$ 

( Samples) &"HZX" 

830  OUTPUT  Vol tmeter ; "Bl, 0B2, "&VAL$ (Points-1) &"X" 

!    Start  and  end  buffer   pointers 

840  OUTPUT  Voltmeter ;"C2G1M2R0X"       !    format, autorange 

850  OUTPUT  Vol tmeter; "NO, "&VAL$ (Points) &"S1,"&VAL$ 

(Samples)  &"HZX" 

86  0  OUTPUT  Vol tmeter ; "Bl ,0B2 , "&VAL$ (Points-1) &"X" 

!  Start  and  end  buffer  pointers 

87  0  PRINT    "WAITING    FOR  BUFFER   TO    FILL    " 
880  OUTPUT  Voltmeter ;"C1T3X" 
890  OUTPUT  Voltmeter ;"C2T25X" 
900  TRIGGER  Voltmeter 
910  WAIT  Delay 
920  OUTPUT  Voltmeter; "C1F0X" 
93  0  PRINT    "CHANNEL   1    READINGS" 
940  FOR  J   =   0   TO   Points-1 
950  ENTER  Voltmeter ;A( J) 

96  0  NEXT  J 

97  0  OUTPUT  Voltmeter; "C2F0X" 
9  80  PRINT    "CHANNEL   2    READINGS" 
990  FOR  J   =   0   TO   Points-1 
1000  ENTER  Voltmeter ;B( J) 
1010  NEXT   J 

1040  OUTPUT  Voltmeter ;"U6X"           ! 

1050  ENTER  Voltmeter ;D$                     ! 

1070  Samples  =  VAL(D$ [3 ,11] ) 

1080  OUTPUT    @Of ile;Points,Freq(I) , Samples, A(*) ,B(*) 

1090  OUTPUT  Voltmeter ; "C1F1XC2F1X"    !    Cancel   waveform 

!   output 

1100  NEXT    I 

1110  ASSIGN   @Ofile   TO   *             !    Close   data   file 

1120  DISP    "DATA    COLLECTION    COMPLETED" 

113  0  END 


Trigger 
Trigger 
Trigger 


source 
source 
both  channels 


Obtain   sampling 
frequency 
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APPENDIX  D 
Source  code  listings 


Dl 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

26  0 

290 

300 

310 

320 

330 

335 

340 

345 

3  50 

370 

3  80 

390 

400 

410 

420 

43  0 

440 

450 

460 

47  0 

480 

490 

500 

510 

520 


********************** 

AUTOMATED   TRANSFER 
This   program  determi 
circuits  and   devices 
lock-in   systems.    The 
the   input  and   referen 
system.    The  Keithley 
output  voltages   of   th 
both   autosensitivity 

The  following  hardwar 
HP8904A   Ch  A: 
HP8904A  Ch  B: 


******************************* 

FUNCTION  MEASUREMENTS  * 
nes  the  frequency  response  of  * 
using  computer  controlled  * 
HP8904A  Synthesizer  provides  * 
ce  signals  of  the  lock-in  * 
194A  measures  the  input  and  * 
e  system.  The  program  provides* 
and  autophase  selection.      * 


e  connections  are  necessary:   * 
Input  signal  * 

Reference  signal  * 


Keithley  194A  Ch  1: 
Keithley  194A  Ch  2: 


Output  voltage 
Input  voltage 


6236B  Power  supply:    +15  V  &  -15  V 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
***************************************************** 


Work  station:   HP9836 
Language:  HPBASIC  3.0 

Programmer:  T.R. Ramesh 

File  name:  ATFM 


03/27/89 


PROGRAM   DECLARATIONS 


OPTION  BASE  0 

INTEGER  Voltmeter,  Fn_ generator,  Printer 

INTEGER  Ref,  Second,  I,J,Z 

DIM  V_in(100)  ,  V_out(100),  Db_  gain  (100)  ,  Freq(lOO) 

DIM   Phase(lOO)  ,A(10)  ,B(10)  ,Offset(100) 

DIM  Mag$[8] ,Phase$[8] ,Name$[20] ,Title$[65] ,Date$[15] 


Voltmeter   =  708 
Fn_generator  =  717 
Printer   =  703 
Second  =  1 
Five   sec   =   5 


HP- IB 
HP- IB 
HP- IB 
Delay 
Delay 


Address 
Address 
Address 


for 
for 
for 


194A 

HP8904A 

HP2631B 


!  Set  up  a  frequency  array  from  10  Hz  to  100  KHz  in   I 
!  logarthmic  steps  ! 

i 1 

•  • 

GOSUB    Initialize  !    Initialize  194A  &  HP8904A 

DATA  10,20,30,40,50,60,70,80,90 

J=0 

FOR  Z=0  TO  4 

FOR  1=0  TO  8 

READ  Freq(I+J) 
Freq(I+J)=Freq(I+J)*10*Z 

D2 


Main   data  acquisition   loop 


530  IF   Freq(I+J)    >   Value   THEN  GOTO  Again 

540  Points=Points+l 

550  NEXT    I 

56  0  J=J+9 

57  0  RESTORE 
580  NEXT   Z 
590 
600 
610 

620  Again:    OUTPUT  Voltmeter ; nC2Xn         !    Point   to   input   of 

630  !    lock-in   system 

6  40  ON   ERROR  GOTO   Errors 

6  50  FOR   I   =    0    TO   Point s-1 

660  DISP   "Data   acquisition   is   occuring   at   this   time" 

670  DISP   "    ("&VAL$(Freq(I) )&"    Hz)" 

680  BEEP 

690  OUTPUT  Fn_ generator ;  "FRA"&VAL$  (Freq (I ) )&"HZ n 

700  OUTPUT   Fn_generator;"FRB"&VAL$(Freq(I)  )  &nHZ" 

710  OUTPUT  Fn_generator;"APA"&Ampa$&"VL  APB "&Ampb$& "VL" 

720  ! 

730  !    Program  Sampling   rate  and  the  number   of   samples   for 

740  !   measurement 

750  IF   Freq(I)<100    THEN    Samples   =   l.E+3 

760  IF   Freq(I)>=100    AND    Freq(I)<lE3    THEN   Samples=l . E+4 

770  IF   Freq(I)>=lE3    AND   Freq(I)<lE4    THEN    Samples=l . E+5 

780  IF   Freq(I)>=lE4    THEN   Samples   =   5. E+5 

790  Num  =   100*Samples/Freq (I) 

800  Here:    OUTPUT   Fn_ generator ; "PHB0DG"       !    Reset   reference 

805  !    phase 

810  OUTPUT  Voltmeter;  "NO,  "&VAL$ (Num) &"XS1,  "&VAL$ 

(Samples)&"X" 

820  OUTPUT  Voltmeter ;"C2F2T3X"         !    RMS,    GET   Trigger 

83  0  TRIGGER  Voltmeter                            !   Trigger 

840  WAIT  Second 

850  ENTER  Voltmeter ;V_ in (I )                !    Input    reading    (RMS) 

860  OUTPUT  Voltmeter ;"ClZOX"              I    Zero   off 

870  OUTPUT  Voltmeter ; "C1XN0 , "&VAL$ (Num) &"XS1 , "&VAL$ ( 

2*Samples)&"X" 

880  OUTPUT   Fn_generator;"APA0VL  APBOVL" 

890  OUTPUT  Voltmeter ;"C1T3X" 

900  TRIGGER  Voltmeter 

910  WAIT  Second 

920  ENTER  Vol tmeter ;Of f set (I )           !    Obtain   offset 

930  OUTPUT  Voltmeter ;"C1Z4X"             !    Enable  zero 

940  GOSUB   Autoset                                       !    Determine   phase 

shift 

950  OUTPUT   Fn_generator;"PHB"&VAL$ (Middle (I) ) &"DG" 

!    Shift   reference   phase  by 
!    the   phase_shift 

980  OUTPUT  Voltmeter ;"C1T3X" 

990  TRIGGER  Voltmeter 
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1000  WAIT  Second 

1010  ENTER  Voltmeter ;V_out (I)   !  Output  reading  (Avg) 

1020  V_out(I)  =  V_out  (I)/(VAL(Ampb$)  ) 

1030  V_out(I)    =  V_out  (I)*Scale    !    Convert  Avg   to  RMS 

1040  Db_gain(I)    =   20*  (LGT(  (  (V_out  (I)  )  )/(V_in  (I)  )  )  ) 

!  Compute  gain 

1050  OUTPUT  Voltmeter ;nC2Xn     !  Point  to  input 

106  0  NEXT  I 

1090   ! ! 

1100  !  Identify  whether  the  phase  angle  is  lead  or  lag    ! 

1110   ! ! 

1120  REDIM  Phase (Points-1) , Offset (Points-1) 

1130  Min  =  MIN(Phase(*) ) 

1140  Max  =  MAX (Phase (*) ) 

1150  Angle$="LAG" 

1170  Off  =  0.0 

1180  FOR  1=0  TO  Points-1 

1190  Off  =  Off+Offset(I) 

1200  NEXT  I 

1210  Off  =  Off/Points 

1220  IF  Min  >  181.0  THEN 

123  0  FOR  I  =0  TO  Points-1 

1240  Phase(I)=Phase(I)-360.0 

1250  Angle $="LEAD" 

126  0  NEXT  I 

1270  END  IF 

1280   ! ! 

1290  !  Obtain  parameters  from  the  user                  ! 

13  00   ! ! 

1310  DISP    "Create   Magnitude  and   Phase   data  files"; 

13  20  INPUT  Ans$ 

1330  Ans$=UPC$(Ans$) 

1340  IF  Ans$="Y"   THEN 

1350  DISP  "Magnitude  data  file  name"; 

136  0  INPUT  Mag$ 

137  0  DISP    "Phase   data  file   name"; 
13  80  INPUT  Phase$ 

13  90  END   IF 

1400  DISP  "Enter  your  name"; 

1410  INPUT  Name$ 

1420  DISP  "Enter  the  title  of  the  experiment"; 

1430  INPUT  Title$ 

1440  DISP   "Enter    today's   date"; 

1450  INPUT   Date$ 

1 46  0      ! ! 

147  0  !  Create  Magnitude  and  Phase  data  files             ! 

1480   ! ! 

1490  REDIM  Freq(Points-l) ,Db_ gain (Points-1) , 

Phase (Points-1) 

1510  IF  Ans$="Y"  THEN 
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Print  data  and  results 


1520  Create_file(Points, Freq(*) ,Db_gain(*) ,Mag$) 

153  0  Create_f ile (Points, Freq (*) , Phase (*) ,  Phase $) 

1540  END   IF 
1550 

156  0 

157  0 

1580  DISP  "Print  data  and  results"; 

1590  INPUT  Print$ 

1600  Print$=UPC$(Print$) 

1610  IF  Print$="Y"  THEN 

1620  GOSUB    Print_header 

1630  FOR    I    =    0    TO   Points-1 

1640  PRINT  USING    " 6D. Df 5X, 3D. 4D, 8X, 3D. 4D, 5X, 3D. 4D, 

4X,3D.4D";Freq(I)  ,V_in(I)  ,V_out(I)  ,Db_gain(I)  , 

Phase (I) 

16  50  NEXT  I 

1660  END  IF 

167  0   GOTO  Fin 

i 


Initialize   Fn_ generator  and  Voltmeter  ! 
i 


1680 

1690 

1700 

1710  Initialize:    ! 

1720  DISP   "Initializing  Voltmeter. . ." 

1730  REMOTE  Voltmeter                 !    Place  194A  in   remote 

1740  CLEAR  Voltmeter                    !   Reset  194A 

17  50  WAIT  Second 

1760  OUTPUT  Voltmeter ; "C1T7XC2T7X"    !    Turn   Ch   1    &    Ch    2    off 

1770  OUTPUT  Voltmeter ;"G1R0X"    !    ASCII   format, autorange 

1780  OUTPUT  Voltmeter ;"Q0X"    !   Disable   reading  buffer 

1790  DISP    "Initializing   Fn_ generator. .. " 

1800  REMOTE   Fn_generator         !    Place   HP8904A   in   remote 

1810  CLEAR   Fn_generator           !   Reset  HP8904A 

1820  DISP   "Signal    amplitude    (Peak)     [    Output   of    test 
circuit    should  not  exceed  11   V    (Peak)]"; 

1825  INPUT   Ampa$ 

1830  DISP   "Reference   Signal    (1-SI,     (Default)    SQ)"; 

1840  INPUT  Ref 

1850  SELECT   Ref 

1860  CASE  1                                 !    Sinusoidal    reference 

1870  Sig$="SI" 

1880  Value  =   100000.0    !    Maximum   frequency 

1890  Scale  =   SQR(2)         !   Scale  factor 

1900  CASE    ELSE 

1910  Sig$="SQ"                    !    Square  wave   reference 

1920  Value  =   50000.0      !    Maximum  frequency 

1930  Scale  =    PI    *    0.5/SQR(2)     !    Scale   factor 

1940  END   SELECT 

1960  Ampb$="l"                                !    Reference   amplitude    (Peak) 

1970  OUTPUT   Fn_ generator ;"PS  GM0   FCIOF   FC20F" 

!    Preset,    Channel    configuration 

!   mode,    Float    controls   off 
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2000 


2040 


2080 
2090 
2100 
2110 
2120 
213  0 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 

2350 
2360 
2370 
23  80 
23  90 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 


OUTPUT  Fn_generator;"FRA10HZ    APAl,&Ampa$&MVL  PHAODG 

WFASI"  !    Ch  A:    Freq   =  10    Hz,    Amp  =   5   V 

!    Phase  =    0   Deg,    Wave  =    Sine 
OUTPUT   Fn_ generator ;nFRB10 HZ    APB"&Ampb$&"VL    PHBODG 

WFBSQ"  !    Ch  B:    Freq    =   10    Hzf    Amp  =   1   V 

!    Phase   =   0   Deg,    Wave  =    Square 
OUTPUT   Fn_generator;">    >"    !    Simulate  NEXT   Key 
RETURN 


Autosensitivity  and  autophase   selection 


Autoset:    ! 
I 

OUTPUT   Fn_generator;"APA"&Ampa$&,,VL  APBl,&Ampb$&"VL" 
OUTPUT   Fn.generator^PHB^ODG" 

FOR  J  =    0   TO  4  !    Output  voltage    in   steps   of   90 

!    degrees   ref.    channel   phase 

OUTPUT  Voltmeter ;"C1T3X" 

TRIGGER  Voltmeter 

WAIT  Second 

ENTER  Voltmeter ;B( J) 

B(J)    =   B(J)/(VAL(Ampb$) ) 

IF  J  =  4  THEN  GOTO  Range 

IF  J  =  3  THEN 

OUTPUT  Fn_ generator, -MPHB3 59. 9DG" 

ELSE 

OUTPUT  Fn_generator;"PHBUP" 

END   IF 
NEXT   J 
t 

Range:    !    Identify   the   range   during  which  quadrature 
phase  occurs 
j 

Equal    =    0 

FOR  J   =   0    TO  3 

IF    SGN(B(J) )=SGN(B(J+1) )    THEN 

Equal  =   Equal   +  1 
ELSE 

IF   SGN(B(J)  )<0.0    THEN 
First   =   J*90.0 
GOTO   Find 
END    IF 
END    IF 
NEXT   J 

IF  Equal    =    4    THEN  GOTO   Senst 
Find:    ! 
Last=First+90.0 
IF   Last=360.0    THEN  Last=359.9 
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2520  !  Find  quadrature  phase  using  Bisection  Search  method 

2  53  0  ! 

2540  OUTPUT  Fn_generator ; "PHB "&VAL$ (Fi rst) &nDGn 

2550  OUTPUT  Vol tmeter ; "C1T3X" 

2560  TRIGGER  Voltmeter 

2570  WAIT   Second 

2580  ENTER  Vol tmeter ;A (0) 

2590  A(0)=A(0)/VAL(Ampb$) ) 

2600  OUTPUT   Fn_generator;nPHBn&VAL$(Last)&nDG" 

2610  OUTPUT  Voltmeter ;"C1T3XH 

2620  TRIGGER  Voltmeter 

2630  WAIT  Second 

2640  ENTER  Voltmeter ;A ( 2) 

2650  A(2)=A(2)/(VAL(Ampb$) ) 

2660  IF  (A(0)*A(2) )>0.0)  THEN  GOTO  Here 

2670  IF  (A(0)<=0.0)  THEN 

2680  Rtbis=First 

2690  Dx=Last-First 

2700  ELSE 

2710  Rtbis=Last 

2720  Dx=First-Last 

2730  END   IF 

2740  FOR  J=l    TO   40 

2750  Dx=Dx*0.5 

2760  Xmid=Rtbis+Dx 

2770  OUTPUT  Fn_ generator ; "PHB "&VAL$ (Xmid) &nDG" 

27  80  OUTPUT  Vol tmeter ; "C1T3X" 

2790  TRIGGER  Voltmeter 

2  800  WAIT  Second 

2  810  ENTER  Voltmeter ;A (1) 

2820  A(l)=A(l)/(VAL(Ampb$)) 

2830  IF    (A(1)<=0.0)    THEN   Rtbis=Xmid 

2840  IF    (ABS(Dx)<0.5)    THEN  GOTO    Found 

2  850  NEXT   J 

2860  Found:    !    Obtain   the   phase_shift 

2870  Middle (I)    =   Rtbis+90.0 

2890  IF  Middle(I)>359.9    THEN   Middle (I) =Middle (I ) -360.0 

2  900  GOTO   Al 

2910  Senst:    !    Sensitivity   variations 

2930  IF  VAL(Ampb$)<10    THEN 

2940  Ampb$    =  VAL$(VAL (Ampb$) +1) 

2950  ELSE 

2960  IF  VAL(Ampa$)<10    THEN 

2970  Ampa$=VAL$(VAL(Ampa$)  +1) 

2980  ELSE 

2990  DISP  "Unable  to  detect  phase  ! ! ! ■ 

3000  STOP 

3010  END  IF 

3020  END  IF 

3  03  0  GOTO  Here 
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3040 

All    RETURN 

3050 

t 

3060 

Errors:    ! 

3070 

IF    ERRN=28    THEN 

3080 

DISP    "NEGATIVE  VALUE    FOUND    !  !  !  ■ 

3090 

OUTPUT   Fn_generator;"PHBODG" 

3100 

GOTO  Here 

3110 

ELSE 

3120 

DISP   "Error"   ERRN    "has  occured"; 

3130 

STOP 

3140 

END    IF 

3150 

GOTO    Fin 

3160 

Print_ header:    ! 

317  0 

I 

3180 

DISP   "Select   a   printing  device    (1-CRT, 2-PRINTER) " ; 

3190 

INPUT  Device 

3200 

IF  Device=2   THEN 

3210 

PRINTER   IS    Printer 

3220 

ELSE 

3  23  0 

PRINTER    IS    CRT 

3240 
3250 

END    IF 

3  26  0 

PRINT   Title$ 

3  27  0 

PRINT 

3280 
3290 

PRINT  Name$ 

3300 

PRINT    "Frequency   range:    10    Hz   TO  100   KHz      Data   points: 

"; Points 

3310 

PRINT   "Output   phase   angle:    ",Angle$ 

3320 

IF  Ans$="Y"    THEN 

3330 

PRINT    "Magnitude   data   file:    ",Mag$,"      Off set:", Off 

3340 

PRINT    "Phase   data  file           :    ",Phase$ 

3350 
3360 

END   IF 

3370 

PRINT 

3380 

PRINT   "Frequency      Input  voltage     Output  voltage 

Db_gain   Phase    " 

3390 

PRINT    "       (Hz)                          (RMS)                            (RMS) 

(dB)         (Deg)" 

3400 

PRINT 

3410 

RETURN 

3420 

I 

3430 

Fin:    PRINTER    IS    CRT 

3440 

END 

3450 

i 

3  46  0 

SUB   Create_file(Num,X(*) ,Y(*) ,Plot_title$) 

3  47  0 

!   Num:      Number    of    data   points 

3480 

!    X(*) :    X  axis   data 

3490 

!    Y(*) :    y   axis   data 

3500 

!    Plot_title$:    Data   file   name 
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3  510  OPTION   BASE    0 

3  520  INTEGER    I, J 

3530  Records=INT( (Num+1) /16) +1 

3540  CREATE  BDAT   Plot_ti tle$, Records 

3550  ASSIGN   @File  TO   Plot_title$ 

3560  OUTPUT    @File ;Num;X (*) , Y (* ) 

3  57  0  ASSIGN   @File   TO    * 

3  5  80  SUBEND 
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Table  Dl 

Sample  output  listing  of  transfer  function  measurement 

using  lock-in  systems 


FREQUENCY 

RESPONSE   OF   A   RC 

LOWPASS    FILTER 

T.R.RAMESH 

25   March  1989 

Frequency 

range:    10    Hz    to   100   KHz 

Data   points:    37 

Output   phase   angle:    LEAD 

Offset: 

-0.0372  V 

Magnitude 

data   file:    LPMAG 

Phase   data 

file           :    LPPHA 

Frequency 

Input  voltage 

Output  volta 

ge     Db  gain 

Phase 

(Hz) 

(RMS) 

(RMS) 

(dB) 

(Deg) 

10.0 

3.5320 

3.2691 

-.6718 

-1.4063 

20.0 

3.5320 

3.2617 

-.6916 

-2.9004 

30.0 

3.5320 

3.268  9 

-.6724 

-4.3006 

40.0 

3.5320 

3.2601 

-.6957 

-5.8008 

50.0 

3.53  20 

3.2495 

-.7242 

-7.2070 

60.0 

3.5310 

3.2370 

-.7551 

-8.6133 

70.0 

3.5310 

3.2226 

-.7939 

-10.0195 

80.0 

3.5310 

3.2065 

-.8374 

-11.4258 

90.0 

3.5310 

3.1887 

-.8857 

-12.8320 

100.0 

3.5310 

3.1719 

-.9315 

-14.1064 

200.0 

3.5300 

2.9147 

-1.6635 

-26.6309 

300.0 

3.5290 

2.6024 

-2.6457 

-36.8262 

400.0 

3.5280 

2.2982 

-3.7228 

-44.8682 

500.0 

3.5270 

2.0296 

-4.7999 

-51.1084 

600.0 

3.5270 

1.8021 

-5.8326 

-56.03  03 

700.0 

3.5250 

1.6119 

-6.7964 

-59.8096 

800.0 

3.5250 

1.4529 

-7.6983 

-63.0176 

900.0 

3.5250 

1.3192 

-8 .5371 

-65.6104 

1000.0 

3.5250 

1.2171 

-9.2367 

-67.8  076 

2000.0 

3.5240 

.6482 

-14.7063 

-78.7  93  9 

3000.0 

3.5250 

.4487 

-17.9049 

-81.5625 

4000.0 

3.5250 

.3073 

-21.1915 

-83.3203 

5000.0 

3.5260 

.2434 

-23.2202 

-8  4.418  9 

6000.0 

3.5260 

.2059 

-24.6705 

-8  4.7  266 

7000.0 

3.5280 

.1759 

-26.0475 

-85.1221 

8000.0 

3.5290 

.1542 

-27.1937 

-85.3418 

9000.0 

3.5300 

.1382 

-28  .1429 

-85.4736 

10000.0 

3.5440 

.1247 

-29.0735 

-85.5176 

20000.0 

3.5490 

.0618 

-35.1825 

-85.0342 

30000.0 

3.5330 

.0411 

-38.6831 

-84.4189 

40000.0 

3.5200 

.0312 

-41.0340 

-81.7383 

50000.0 

3.5060 

.0256 

-42.7370 

-81.5186 
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ft  ***********************************  ******  icmm  r>  *  t  ttttt  ******* 
* 

*  SOURCE  FILE:     ;:fer_fn.c 
* 

* 

*  FUNCTION:        7.3in() 

* 
* 
* 

*  DESCRIPTION:     Obtains  the  freauency  response  of  s  linear 

*  system  by  the  following  methods. 

*  1.  Three  parameter  lesst  sous  red  fit 

*  algorithm 

*  2.  Cross-correlation  method 
* 

* 

*  DOCUMENTATION 

*  FILES:  atfm.d 
* 

* 

*  ARGUMENTS:       None 
* 

* 

*  RETURN:  None 
* 

* 

*  FUNCTIONS 

*  CALLED:         corr()»fU()»lsf() 
* 

* 

*  AUTHOR:  T.R.Ramesh 
* 

* 

*  DATE  CREATED:     26  March  1999  Version  1.00 
* 

* 
************************************** ******* ttttttttttttttttttttf 

♦include  'atfm.h'  /*  Contains  error-codes  and     %t 

ft   definitions  used  by  xfsr_?n  •* 

/*  routines  \f 


/*  General  purpose  indices      %•' 

ft    Number  of  frequencies 

/*  Selection  variable  I 

/*  Index  for  freauency  I  cop 

/*  Points  for  esch  frwiauency    1/ 

double     sainCMAX].'  ft   Amplitude  response  anvrj     t/ 

phase.shif tCMAX3i    /*  Phase  response  arr3s         I 

Dll 


main( ) 

■C 

int 

i » c  >  J  >  k  > 

1  e  n  2 1  h  t 

type  i 

f  reQf 

points? 

resolutionCMAX] »     /%  Resolution  in  decrees  $/' 

f reauenciesCMAX] >    /%  Freouency  3rray  %/ 

yinCMAX]?  youtCMAX]>  /*  Input  3nd  output  samples  %/ 

samp_frea»  /*  Sampling  freauency  */ 

xCHAX3i  /*  Time  array  */ 

step*  /%  Step  length  */ 

in_amp»  out_amp»     /%  Amplitude  estimates  %/ 

in.offr  out.offj     /*  Offset  estimates  %/ 

in_phase>  out_phase»  /%  Phase  estimates  */ 

max*  omax?  imaxJ       /#  Dummy  variables  "'■' 

char       dsta_fileC20Jf       /*  Input  data  file  */ 

ma3_fileC203»        /*  Magnitude  data  file         */ 
phase_fileC20]r      /*  Phase  data  file  */ 

COMPLEX    inCMAX3»  /*  Input  signal  array  %/ 

outCMAXD»  /*  Output  signal  sppsa 

resultCMAX]>         /%   Correlation  result  vector    */ 

FILE       * input*  /*  Data  file  pointers  %/ 

*mag» 
*phai 

extern  corr( ) »f ft( ) »lsf < ) i 

/t %i 

/%     Obtain  the  input  and  output  samples  of  the  transfer  function*/ 

/%     measurement  %/ 

f% t/ 

putsC\n  AUTOMATED  TRANSFER  FUNCTION  MEASUREMENTS* )  > 
puts('\n\ri  Refer  to  documentation  file  for  datd  file  format")? 

printf('\n  Select  one  of  the  following  methods* ") 5 
ppintf("\n  It  least  sauared  method')? 
printf('\n   2.  Cross-correlation* ) ? 
printf('\n\n  your  choice?' > J 
scanf<"Zda>  Stspe). 

printf('\n  Data  file  name7')? 
scanf  ("is"-  Sdata.file)? 

input  =  f open( data_f i le i ' r' ) ? 

fscanf (inputi "%d" >  Slength) » 

length-1 5 

for  (frea=0i  free  <  length*  freo+T) 

■C 

fscanf ( input; *%d" >  Spoints ) : 
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fsc3nf  (inputr  '%lf  Zlft    Jf  reauencie-sCf  ree]i  SsaftP.frea  )  : 
for  (i  =  Of  1  <  points?  iff) 

fscanf (input »  "/.Iff    8«inCi3)  5 
for  ( i  =  0  ?  l  <  points?  i  +  + ) 

fscanf( input > 'Zlf ' »  SyoutEij) • 
step  =  1 ,0/53mp_f rea? 
/* %f 

/%      Obtain  amplitude  and  phase  responses  for  each  frequency     %/ 
/% !  , 

switch  (type) 

case   LSFJ  /*   Least    sauared   method  */ 

for    (i    =   0?    i    <   points?    i++) 

xCi]    =    i    *   step? 
if    (sine_f it(points>    «»    sim    Sin_aijiP*    f  reouenciesC  f  re<*j  > 

8in_offi    sin_ph3se)    !=   NORMAL) 

< 

ppintf('\n  Error  in  sine  fit  to  input  data')? 

exit(0)  ? 

if  (sine_f  it(points>  :;>  yout>  &out_anp>  f  r  sauenc iesCf  rexl 
Sout_off>  Sout.Phase)  !=  NORMAL) 
■C 

printf('\n  Error  in  sine  fit  to  output  data*)? 
e::it(0)  ? 
> 
sainCfrea]  =  20  *  Io2l0(out_3mp/iri_ai?ip  )  ? 

if  ((out-Phase  -  in.phase)  >  180.0) 

in_phase  +=  360.0? 
else 

■C 

if  ((in.phase  -  out.phsse)  >  130-0) 
out.phsse  +=  360.0? 

■\ 

j" 

phase.shif tCf rec<]  =  out.phsse  -  in_pha;e? 
break? 


case  C0RR: 
points  =3? 
resolutionCf reo]  =  360 . 0/( double )poin ts? 


/%   Cross  correlation  aethod   K/ 


for  (i  s  Of  i  <  points?  i++) 
{ 

inCi]  =  chipIx (yinCi ] > ZERO)  ? 
OutCi]  =  cnplx  (youtCi]  >  ZERO)  ? 
> 
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if  (corr(in>out>  points*  result)  !=  NORMAL) 

■C 

priritf('\n  Error   in  correlation    routine')* 

exit CO) i 

> 
max  =  0.0  J 
for  (i  =  0?  i  <  points?  i++) 

if  (  resultCi]  .  re  >  max)  max   =  resultCi] . re? 

i  =  Oi 

while  (i  <  points) 

■C 

if  (max  ==  resultCi] . re)  break? 

i  +=  it 

y 

phase.shiftCf reo]  =  (double)i  *  resolutior.Cf  reel  ? 

if  (i  >=  points/2) 

phase_shiftCf  rec<]   =  phase_shiftCf pea]   -  160.0' 

phase.shiftCfrea]   *=   -1,0? 

imax   =  0.0? 

0JII3X    =    0.0? 

for  (i  =  0?  i  <  points?  i++) 
■C 

if  (yinCi]  >  imax)  imax  =  yinCi]? 
if  (youtCi]  >  oir.ax)  oinax  =  youtCi]? 

SainCfreo]  =  20.0  *  los:10(  om3;:/iin .?;■;)  ? 
break? 


default! 

exit (0)5 


/*  Create  ssin  and  phase  data  files 


prmtfCXn  Create  data  files?  (Enter  1)')? 
scanf  ( 'Xd'rSc)? 

if  (c  ==  1) 

J" 
!■ 

printf('\n   Gain   data   file?')? 
scanf <,Zs">    Smas.file) ? 
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ppintf(*\n  Phase  dete  file?,)> 
scanf  <  '7.s'  r    8phase«fileJ  ? 
mag  =  f open ( maS-f 1 le« * w " ) » 
phs  =  f openi ph<3se_f  i ie t  ' w"  5  I 
fprintf (ms3> 'Zd' >  length) ? 
fprintf  (pilar' Zdai  length)  J 

for  (i  =  Oi  i  •■'  length}  i++l 
{ 

fprintf  (rnsSf  '\n%f  7.T    '  t  f  recuenciesE  1 3  >  3ainCi3)i 
fp rintf (Ph3> 'NnZf  Xf p f reuuenciesCi] rPhase.shif tCi] ) J 
> 


D15 


*  SOURCE  FILE:     corr.c 

*  function: 
* 

*  description: 

* 

*  documentation 

*  files: 

* 

*  arguments: 

* 

v_in 


corr(v_in>  v_out>  n»  result) 

Determines  the  correlation  of  two  complex 
arrays 

None 


* 

* 
* 
* 
* 
* 
* 
* 

* 
* 
* 
* 


v_out 


result 


(input)  COMPLEX 

Array  containing  the  input  signal 


(input)  COMPLEX 

Array  containing  the  output  signal 


(input)  int 

Size  of  the  array,  (must  be  3  -ower  of  2) 


(output)  COMPLEX 

Array  containing  the  correlation  of  the 

two  input  arrays 


*  return: 

* 

*  functions 

*  called: 

* 
* 

*  author: 
* 

*   DATE  CREATED: 

* 


int 

NORMAL     :  Normal  return 

ERR.INDFT   !  Error  in  input  voltage  DFT 

ERR-OUTDFT  '.  Error  in  output  voltage  DFT 

ERR.PRIDFT  !  Error  in  product  IDFT 


sin( ) »cos( ) 

caddOf  csub()»  cmult'.  )♦ 


T . R.Ramesh 

23  Mar  193? 
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Version  1.00 


*   REVISIONS:       None 
* 

♦  include  'atfm.h'  /*  Header  file  for  ;:fer_fn      %/ 

/%    routines  1/ 

int  corr(v_irif  v_out>  n»  result) 

int        nJ 

COMPLEX    v_inC]» 
v_outC] t 
resulted 

{ 

int        i?Jf  (%   General  purpose  indices      %/ 

extern   int  fft( )  i 

/% 1/ 

/%     Compute  the  DFT  of  both  the  input  vectors  %/ 

/% %/ 

if  (fft(V-in»  n»  DFT)  !=  NORMAL) 
return  (ERR_INDFT)J 

if  (fft(v_OUti  r.F  DFT)  !=  NORMAL) 

return  (ERR.OUTDFT) i 
/% */ 

/%     Obtain  the  product  of  one  vector  with  the  complex         %/ 

/%     conjugate  of  the  other  */ 

/% x/ 

for  (i  =  0  5  i  <  n>  i++) 
■C 

v.inCi]  =  cconJ( v_inCi] ) \ 
resultCi]  =  cmult ( v_inCi]> v.outCi J) i 

/% 1/ 

/%     Obtain  the  IDFT  of  the  resultant  vector 

/% %/ 

if  (fft( result.  n>  IDFT)  !=  NORMAL) 

return  (ERR.PRIDFT)  ; 
/* $/ 

ft      Normal  termination  " ' 

1% */ 

return  (NORMAL); 
> 
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* 

*   SOURCE  FILE!     dft.c 

* 
* 
*  function: 


fft(dst3»  size>  sign) 


*  description: 
* 

* 
* 

*  documentation 

*  files: 

* 

*  arguments: 

t  dat3 

* 

* 

* 

*  size 
* 

* 

*  si  an 
* 

* 
* 

*  return: 
* 

* 

*  functions 

*  called: 

* 

* 

*  author: 

*  date  created: 

* 

*  revisions: 


Determines  the  Discrete  Fourier  Transform 
(DFT)  or  Inverse  Discrete  Fourier  Transform 
(IDFT)  of  3  complex  array. 


None 


(input/output)  COMPLEX 

Input  data  array,  data  contains  DFT  or  IDFT 

upon  return. 


(input)  int 

Size  of  the  array,  (must  be  a  power  of  2) 


(input)  int 

DFT   :   Discrete  Fourier  Transform. 

IDFT  :   Inverse  Discrete  Fourier  Transform. 


int 

NORMAL   :  Normal  return 
ERR.FFT  :  Number  of  data  points  is  not  3 
power  of  2, 


sin( ) «cos( ) 

cadd()>    csubOi    cnultO* 


T  .R.Ramesh 


20    Msr    1989 


None 


Version   1 .00 


D18 


♦include  'atfm.h' 


/*  Header  file  for  xfer.fn      */ 

/*  routines  %.' 


int  fft(data»  size>  sian) 

int      size> 
sian* 

COMPLEX  dstsHi 


■C 
int 


i  f  j  j  k  r  1 1 

mm  a;:  > 

step? 


COMPLEX   wi 

dummy  *  temp? 


/*  General  purpose  indices  */ 

/*  Array  size  for  recursive  # 

/*  transform  computations  */ 

ft   Step  size  to  set-ess  groups  */ 

/*  of  arrays  */ 


/*  A  comple>:  number 
/*  Dummy  variables 


*/ 
*/ 


/%      Check  if  the  size  of  the  data  array  is  a  power  of  2  X/ 

/% */ 

if  ((size  V.      2)  !=  NORMAL) 

return (ERR_FFT) > 
k  =  size? 


while  (k  >1) 
•C 
if  (<k  V.    2)  !=  NORMAL) 

return  (ERR.FFT)i 
k  =  k/2» 
> 


/*  Perform  bit  reversal  of  the  input  data  stream  */ 

f% » 

j  =  0> 

for    (i    =   Of    i    <   size?    i++) 
•C 
if    (l    <    J) 

i  /%   Swap   the    two   conple."  i/ 

dummy. re    =    dstaCiJ.re?  /*   numbers  I 

dummy. lm    =    dataCiD.imJ 
dataCil.re    =    dataCJD.re? 
d a t a C i  3  .  i 7i    =   dstsCJD.iiaJ 
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data C J]. re  =  dummy. re? 
d s t a C  J 1 .  i hi  =  dummy . im? 

> 

k  =  size/25 

while  (  k  >=  1  SS  k  •=  J) 

■C 

J  -=  kJ 

k  =  k/2J 

> 
J  +=  k? 


/*- 
/% 

/* 
/*• 


Perform  FFT  usina  the  Danielson-Lanczos  algorithm 
(successive  doubling  method) 

mma>:  =  1 ? 

while  (size  >=  (2  *  mmsx) ) 
■C 

step  =  2  *  mmex? 
dummy  =  cmp la < ONE t ZERO >» 
w  =  cmpl>:(cos(   PI/(  (douhle)mms:<  ))»  sin(  PI/ 

( (double) (sian  t   mmax  ))))? 
for  (J  =  0?  J  <  mm3>:?  J  f+) 
■C 

for  (1  =  J?  1  <  size?  1  +=  step) 
■C 

temp  =  cnult (dstsCl+mmsx] >  dummy)? 
dataCl-f-mma::]  =  csub(dstsCl] .»  teoiP>? 
d3teCl]  =  C3dd(d3t3Cl]f  tern?)? 
> 
dummy  =  cmult(w»  dummy)? 
j 
mmax  =  step? 


•*/ 
*/ 
*/ 

-  * .  ■ 


/* %/ 

ft     Divide  the  DFT  by  the  size  of  the  array  to  obtain  IDFT      */ 
/* %/ 

if  (si2n  ==  IDFT) 
< 
for(i  =  0?  l  <  size?  i++) 

dated!)  =  cdiv(dstaCi]  •  cmpl::(  (double)size>  ZERO))? 
> 

./* */ 

/%      Normal  termination  */ 

./* */ 

return(NORMAL) ? 

> 
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* 

*  SOURCE   file:  lsf.C 
* 

* 

*  FUNCTION:       sine_f it(np»  >c_dat>  a-dsti  saplitudei  ?-ia> 

*  dc»  thets) 
* 

* 

I  DESCRIPTION:     Provides  the  amplitude*  phase  and  offset 

*  estimates  of  3  sampled  sinewave  us  ins  the 

*  three  parameter  (Known  freauency)  least 

*  sauared  method. 
* 

* 

*  DOCUMENTATION 

*  FILES:  None 
* 

*  arguments: 
* 

*  np  (input)   int 

*  Number  of  data  points. 
* 

*  >;_dat         (input)   double  * 

*  Abscissa  values  of  the  sampled  sisnal. 
* 

t 

*  y_dat         (input)   double  * 

*  Ordinate  values  of  the  sampled  signal. 

*  amplitude     (output)   double  J 

*  Amplitude  estimate  of  the  sampled  signal. 
* 

* 

*  frea         (input)   double 

*  Freauency  (known)  of  the  sampled  signal. 

* 

*  dc  (output)   double  * 

*  dc  offset  estimate  of  the  sampled  signal. 

* 

*  theta         (output)   double  t 

*  Phase  ansle  estimate  of  the  sampled  signal. 

*  FUNCTIONS 

*  CALLED:  cosO?    sin<) 
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*  AUTHOR:  T.R.Ramesh 

* 

*  DATE  CREATED:    21  Feb  1939  Version  1.00 
* 

**********************************  *************  *****  *************/ 

♦  include  'atfm.h' 


int   sine_fit(np*    x_datr    y_dst>    amplitude.'    freaj    do    theta) 

/*  Number  of  points  sampled  */ 

/*  Abscissa  values  */ 

/*  Ordinate  values  */ 

/*  Amplitude  estimate  t/ 

/*  Freauency  )'■/ 

/*  Offset  estimate  */ 

/*  Phase  ansle  estimate  */ 


int 

np  r 

double 

K_datC]» 

y_datC3» 

*amplituder 

f  reaj 

*dc> 

*theta  i 

int     i?Jrk>  /*  General  purpose  indicies  */ 

/*  Dummy  variables  */' 

double  max*  min»  amp_est»  off.estf  sisri*  analer  phif  w* 

sum_yn»  sum_an»  sum_bn>  sum_abri>  sum_aaru  sum-bon* 
sum.yarif  sum_ybn»  sum_yyru  alpha*  betar  vbar*  alpha-bar? 
beta_bar>  a_n>  a_d»  b_n?  b_d?  a»  b»  c5 

/* %/ 

/*  Find  the  initial  amplitude!  offset  and  phase  estimates     */ 
/% */ 

min  =  100.05  /*  Determine  maximum  and    */ 

max  =  0.05  /*  minimum  values  in  the    %/ 

/*  sampled  sisinsl  array     */ 
for  (i  =  Or    i  <  np»  i++) 

< 

if  (y.datCi]  >  max)  max  =  y_datti]5 

if  (y.datCi]  <  min)  min  =  y_dat£i]5 

> 
amp_est  =  (max  -  min)  /  2.0  5 
off.est  =  (max  +  min)  /  2.0? 

w  =  2.0  *  PI  *  freai 
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/% 1 

/*   Estimate  parameters  usina  the  three  psrsmeter  lesst        */ 
/*  souared  fit  technioue  1/ 

/* ■■:/ 

sum_yn  =  ZERO?  /"*  Initialize  sums        */ 

sum_an  =  ZERO? 
sum_bn  =  ZERO? 

sum_abn  =  ZERO? 
sum.aan  =  ZERO* 
sum_bbn  =  ZERO? 

sum.aan  =  ZERO* 
sum.ybn  =  ZERO' 
sum_yyn  =  ZERO? 

ft      Compute  nine  sums  reauired  for  the  estimates  */ 

for  (i  =  0?  i  <  np»  iff' 
-C 

alpha  =  cos<w  *  K.dat£i3)f 
beta  =  sin(w  *  x„d3tEi])5 

sum_yn  =  sum.yn  +  y_datCil» 
sijm.an  =  sum_an  +  alpha ? 
sum_bn  =  sun.bn  +  beta? 

sum.abn  =  sum_abn  +  3 If  ha  *  beta? 
sum_aan  =  sum.aan  f  alpha  *  il?ns? 
sum_bbn  =  sum_bbn  +  beta  %   beta' 

sum_yan  =  sum_yan  +  y.detCiZ  %    alpha? 
sum_ybn  =  sum_ybn  +  y_ds*Ci]  X  b=rtar 

sum.yyn  =  sum_yyn  +  y_datCi3  %    a_d<:tCi3  5 

■l 

/%      Compute   the   following   parameters   using    the    »uas  *' 

/%      calculated   3bove  '-■  ■' 

ybar    =    (sum_yn    /    ( (double)np) )  • 
3lpha_bar   =    (sum_3n   /    ( (double)np) ) 5 

bet3_bar   =    (sum_bn   /    ( (double )np) ) ? 

3_n  =  (((sum_yan  -  ybar  *  sum_an)  /  (sum_sbn  -  betd_b«r  * 
sum. an))  -  ((sum_ybn  -  -ihs?  %  sum_bn)  /  (sun.bbn 
-   beta_bar   *   sum_bn)))? 

s_d   =    (((sum.aan   -   alpha. bar   t   sum_an)    /    fsuft.abn   - 

bet3_b3r   *    sum_an)>    -    ( ( sum_abn    -   alpha. bar    * 
sua.bn)    /    (sum_bbn   -   bets-bar    '<   5iim_bn ) ) )  > 
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b_n  =  ( ( ( sijm_a3ri  -  ybsr  *  suiri_sn)  /  (sum_33n  -  3lph3_b3r 

*  sum_3n))  -  ((sum_ybn  -  yb3r  t   5Ui»_bn)  /  (sum-sbn 
-  slphs-bsr  *  sunt-bn) ) )  > 

b_d  =  ( ( (sum_3bri  -  bets.bsr  *  sum.sn)  /  <sum_53n  - 

slph3_b3r  *  susi-sn))  -  ((sum.bbn  -  bets_b3r  * 
sum_bn)  /  (sum-sbn  -  slphs-bsr  H  sum-bn)))? 

3  =  3_n  /  3_d? 
b  =  b_n  /  b.di 

c  =  (yb3r  -  (3  *  3lphs_b3r)  -  (b  *  bets.bsr))* 
Samplitude  =  sart(3  *  3  +  b  *  b)J 
*thet3  =  3tsn(  1.0  *  s/b)> 
*thet3  =  *thet3  *  130.0  /  PI  J 
if  (b  <  0.0)  Kthets  -   Uhets  +  130.0? 
if  (b  >  0.0  II   3   <   0.0)  *theta  =  *thets  +  360.0- 
*dc  =  c> 
/* */ 

/*   Noruisl  terminstion  */ 

f% %/ 

return  (NORMAL) J 
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* 

3tf m .h 


* 
* 

* 

* 
* 
* 
* 

* 
* 

* 


source  file: 


description: 


documentation 
files: 


author: 


Header  files  for  transfer  function 
measurements  using  cross  correlation  and 
three  parameter  least  sauared  fit  algorithm 


None 


T.R.Ramesh 


Version  1 .00 


*   DATE  CREATED:    22  March  198? 

* 

* 


♦include  <stdio.h> 
♦include  <math.h> 
♦  include  *  com:*  lex.  h* 


/*  Header  file  for  complex 
/*  functions 


♦define 

MAX 

1000 

♦define 

PI 

3.141592653 

♦define 

DFT 

1 

♦define 

IDFT 

-1 

♦define 

LSF 

1 

♦define 

CORR 

2 

♦define 

ONE 

1.0 

♦define 

ZERO 

0.0 

/*  Maximum  number  of  points  *.' 

/$  Constant  PI  */ 

/*  Discrete  Fourier  Transform  *, 
/*  Inverse  Discrete  Fourier 
•'*  Transform 

/*  Least  souared  method  */ 

/%   Cross  correlation  method  #/ 


t% %/ 

/%     Error  codes  for  functions  used  in  ::fer_fn  routine  %/ 

/% . 

/*  Normal  return  */ 

/%  Return  error  code  for  ffti)\/ 
/*  Error  in  input  voltoSe  DFT  t/ 
/  *  Error  in  output  vo 1 1 afle  DFT t < 
/*  Error  in  product  IDFT      I 


♦define 

NORMAL 

0 

♦define 

ERR-FFT 

1 

♦define 

ERR_INDFT 

■-> 

♦define 

ERR.0UTDFT 

3 

♦define 

ERR_PRIDFT 

4 
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* 

*  source  file:    atfm.d 

* 

*  description: 
* 
* 

*  documentation 

*  files: 

*  author: 
* 

*   DATE  CREATED: 


Provides  details  regarding  file  formats 
for  :cfer_fn  routines 


None 


T.R.Ramesh 


22  April  198? 


Version  1  *  00 


Transfer  function  measurements  of  linear  satems  are  accomplished 
us  ins2  two  algorithms: 

It  Three  parameter  least  souered  fit  algorithm 
2.  Cross  correlation  of  input  and  output  signals 

The  input  data  file  for-  both  the  methods  fftust  be  of  the 
following  format: 


(i)    Number  of  freauencies  '.  integer) 
For  every  freauencs* 
(l)    Number  of  joints 


( ii )  Freauency 

(iii)  Sampling  freauency 

( iv)  Input  samples 

(v)  Output  samples 


( integer) 
(must  be  an  integer  Fewer-  of  2 
for  cross  correlation  method) 
(double! 
i double 
(double 
(double 


The  output  files  created  sre   of  the  following  format! 
(i)    Number  of  points       (integer; 
For  every  freauency 
(ii)   Freauency  (double'1'  Gain/Phase   (double) 
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ABSTRACT 

This  thesis  discusses  hardware  and  software  techniques 
for  automating  transfer  function  measurements  of  linear 
sytems. 

The  hardware  approach  is  using  computer  controlled  lock- 
in  systems.  The  lock-in  system  is  designed  using  a 
precision  IC  Multiplier  AD534.  Software  techniques 
considered  are  the  three  parameter  (known  frequency)  least 
squares  fit  algorithm  and  the  cross  correlation  of  input 
and  output   signals   of   the   linear    system 

The  hardware  approach  provides  phase  accuracy  dependent 
on  the  resolution  of  the  phase  shifter.  The  amplitude 
response  is  a  function  of  the  phase  response  due  to  the 
nature   of   the  measurement   process. 

Software  techniques  thoroughly  rely  on  the  simultaneous 
sampling  of  input  and  output  signals.  The  three  parameter 
method  is  faster  and  provides  better  accuracy  than  the 
cross  correlation  method  for  a  given  number  of  samples  per 
cycle. 


